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Abstract 



We consider finite temperature correlation functions in massive integrable Quantum Field 
Theory. Using a regularization by putting the system in finite volume, we develop a novel 
approach (based on multi-dimensional residues) to the form factor expansion for thermal 
correlators. The first few terms are obtained explicitly in theories with diagonal scattering. 
We also discuss the validity of the LeClair-Mussardo proposal. 



1 Introduction 



One of the central tasks in many-body quantum physics is the calculation of correlation func- 
tions. They yield considerable amount of information about the dynamics of the system and their 
Fourier modes can be measured, for example with elastic neutron scattering experiments. In ad- 
dition to the correlations in the ground state it is also important to calculate correlation functions 
at finite temperatures, in which case the system is populated by a number of excited states. 

In this paper we consider finite temperature correlation functions in 1+1 dimensional integrable 
(T) ■ models. A common property of these theories is that their Hamiltonian can be diagonalized 

OO . analytically using the Bethe Ansatz [1, 2]. Moreover, there are powerful methods available to 

obtain correlation functions, the most general being the so-called form factor expansion. The 
idea is to obtain the matrix elements of local operators (form factors) in the eigenstate basis of 
the Hamiltonian, and then to sum up the spectral series. The difference as compared to usual 
approaches (like perturbation theory) is that in integrable models both the spectrum and the form 
factors can be calculated exactly. This presents a unique opportunity to study strongly correlated 
quantum systems in situations where conventional methods break down. Interest in integrable 
models has recently been renewed, in large part due to the developments in recent years that 
made it possible to realize certain models with the help of optical and magnetic traps [3, 4, 5, 6]. 

The ideas of the form factor expansion are quite general, however the methods to obtain the 
form factors and to sum up the spectral series can be different. One framework is provided by the 
Algebraic Bethe Ansatz (ABA) [2, 7, 8], which was applied successfully to a number of models, 
most prominently the Heisenberg spin chains [9, 10. 11] and the ID Bose gas [12, 13]. Correlation 
functions are obtained typically in the form of integral series [14, 15] or multiple integral formulas 
[10, 11], or one can resort to numerical summation schemes [16, 17]. Moreover, there are methods 
to handle the finite temperature situation, either through generalizations of the basic techniques [2] 
or developing an alternative description using the so-called Quantum Transfer Matrix [18. 19, 20]. 

Integrable Quantum Field Theory provides a different framework to obtain form factors and 
correlation functions. In these theories the basic object is the factorized S-matrix [21, 22], and 
the relation to a microscopic description is rather indirect. It can be shown, that the form factors 
satisfy a certain set of equations (the form factor bootstrap equations) which follow from general 
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field theoretical arguments supplemented with the special analytic properties of the S-matrix 
[23, 24, 25, 26]. The idea is to provide a general solution to these equations, and then to identify 
those solutions which correspond to a given local operator [27] . The resulting form factor functions 
can then be used to construct correlation functions. An important feature is that the form factors 
are calculated in the infinite volume asymptotic state basis. This is to be contrasted with the 
situation in ABA, where one starts with a finite system and the infinite volume limit is only 
performed afterwards. Interesting connections between the form factor bootstrap and the ABA 
were pointed out recently in [28]. 

The problem of zero-temperature correlations in integrable QFT is well understood. Although 
an analytic summation of the spectral series is in general not possible (except in some simple 
models [29, 30, 31]), the series has very good convergence properties in massive models and can 
be evaluated numerically to a desired precision [25, 32]. On the other hand, the problem of finite 
temperature correlation functions is less understood and it has been subject to an active research 
in the last ten years [33, 34, 35, 36, 37, 38, 39, 40, 41]. The idea is to use the zero-temperature 
form factors in the spectral series, however, the thermal average is ill-defined in infinite volume. 
The problem is related to the appearance of disconnected terms in the expansion, which lead to 
formally divergent expressions. Following Balog it can be shown that the divergent parts cancel 
with contributions from the partition function [42], however it is a highly non-trivial task to obtain 
the finite left-over pieces. 

There have been attempts to write down a regularized version of the spectral series. In 
particular, LeClair and Mussardo proposed an expansion for the one-point and two-point functions 
in terms of form factors dressed by appropriate occupation number factors containing the pseudo- 
energy function from the thermodynamical Bethe Ansatz [34]. Their proposal for the two-point 
function was questioned by Saleur [35] (see also [43]); on the other hand, he also gave a proof of the 
LeClair-Mussardo formula for one-point functions provided the operator considered is the density 
of some local conserved charge. However, it was demonstrated in the case of one-point functions 
that the results obtained by naive regularization are ambiguous; this is shown in particular by the 
difference [38] between the formulae proposed by LeClair and Mussardo and by Delfino [37] 

This motivated the present authors to develop a regularization method based on finite volume 
form factors [44, 45] which was applied to one-point functions giving a confirmation of the LeClair- 
Mussardo formula. The central idea is to use a finite volume setting to regularize the divergences, 
and to compute the physical quantities in finite volume. At the end of the calculation, the 
volume is taken to infinity. If one computes only quantities meaningful in infinite volume, the 
divergences cancel and the end result is well-defined. Because finite volume is not an ad hoc, 
but a physical regulator (note that physically realizable systems are always of finite size) one is 
virtually guaranteed to obtain the correct result provided the calculations are performed correctly. 
The existence of a mass gap m is essential in this approach: the Boltzmann-factor e~ m l T provides 
a natural small parameter for the finite temperature expansion. The result is an integral series, 
where the Ath term represents A-particle processes over the Fock- vacuum. The contributions with 
a low number of particles can be interpreted as disconnected terms of matrix elements calculated 
in a thermal state with a large number of particles [41, 46]. In this sense the approach is similar 
to the one used in the seminal works in ABA [14, 15]. However, a distinctive feature is that 
the calculations do not use any information about the form factors other than their singularity 
properties. 

The finite volume form factor approach was extended to boundary operators as well [47], which 
was used to compute finite temperature one-point functions of boundary operators [48]. Another 
application of the bulk finite volume form factors is the construction of one-point functions of bulk 
operators on a finite interval [49]. One of the present authors have also used this formalism in 
[50] to construct the form factor perturbation expansion in non-integrable field theories (originally 
proposed by Delfino et al. [51]) beyond the lowest order. The finite volume regularization method 
was also exploited in a recent work by Essler and Konik [52, 53]. They computed the first nontrivial 
contribution to the dynamical spin-spin correlation function of a spin chain using the 0(3) sigma 
model. However, the methods they use do not have any obvious extension to higher order, albeit 
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it is reasonably clear that the finite volume regularization must work. 

In this paper we develop a systematic method to compute the finite temperature form factor 
expansion to arbitrary orders, relying on an application of the residue theorem for multiple complex 
variables. We also demonstrate that the same method can be applied to computing the form factor 
perturbation theory contributions, and the zero-temperature three-point function. 

The outline of the paper is as follows: section 2 summarizes the necessary background on form 
factor bootstrap and finite volume form factors. Also, it sets up the framework for the expansion 
of the thermal two-point function, and discusses the scope and validity of the approach. In section 
3 we present the calculation of the (zero-temperature) 3-point functions as a warm-up example. 
The expansion for the finite temperature 2-point function is derived in section 4, except for 
some technical details that are relegated to appendices. The discussion of the LeClair-Mussardo 
proposal in light of our results is presented in subsection 4.5. We present an application of our 
method to form factor perturbation theory in section 5, and conclude in section 6. 



2 Form factor expansion for the thermal two-point function 

2.1 The form factor bootstrap 

Here we give a very brief summary of the equations of the form factor bootstrap, in order to 
set up notations and to provide background for later arguments; the interested reader is referred 
to Smirnov's review [54] for more details. For the sake of simplicity let us suppose that the 
spectrum of the model consists of a single particle mass m. The energy and the momentum of 
an on-shell particle is parametrized by the rapidity variable as E = mcosh0 and p = msinh0. 
Because of integrability, multi-particle scattering amplitudes factorize into the product of pairwise 
two-particle scatterings, which is described by a pure phase, which we denote by S (0) where is 
the relative rapidity of the incoming particles. Incoming and outgoing asymptotic states can be 
distinguished by ordering of the rapidities: 

10 q \_ J |0i>---A) m ■ 9i> 6 2 > ■■■ > 6 n 

a> ou * : e 1 <e 2 <...<e n 

and states which only differ in the order of rapidities are related by 

[0i, . . . , 6k, 9k+i, ... ,9 n ) = S{6k — 6k+i)\6i, . . . , 6k+i,6k, ■ ■ ■ , 9 n ) (2.2) 

from which the S matrix of any multi-particle scattering process can be obtained. The normal- 
ization of these states is specified by the following inner product for the one-particle states: 

(6'\6) = 2tt5(9' - 0) (2.3) 

The form factors of a local operator 0(t,x) are defined as 

10(0,0)10!,... A) (2.4) 

With the help of the crossing relations 

f mn(^l! ■■■■> #ml#li ■ ■ ■ > #n) = F^_ ln+1 (9 1 , . . . ,9 m ^ 1 \9 m + ITT, 01, ... , 0„) 
n k—1 

+ ^5{9' m - 9 k ) J] S(9 t - 9 k )F%_ ln _ 1 (9[,..., 9' m _ x \6 U fe -i, fe+ i ...,6 n ) (2.5) 

k=l 1=1 

all form factors can be expressed in terms of the elementary form factors 

F n °(0!,...,0 n ) = (O|O(O,O)|0!,...,0 n ) (2.6) 

which satisfy the following equations: 
I. Lorentz transformation: 
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F°{d! + A,0 2 + A, . . . ,9 n + A) = exp ( So A) F n °{9 1 ,9 2 , ...,6 n ) (2.7) 

where sq denotes the Lorentz spin of the operator O. 

II. Exchange: 

F n ■ ■ ■ ,dk, @k+l, ■ ■ ■ , On) = 

S(0k — 0k+l)Fn(0l, • • • , Ok+l,&k, ■ ■ ■ , On) (2.8) 

III. Cyclic permutation: 

F%{9 X + 2nr, 9 2 ,...,9 n ) = F°(9 2 , 9 n , 9{) (2.9) 

IV. Kinematical singularity 

-iResF° +2 (9 + m,9',9 1 ,...,9 n ) = ( 1 - TT 5(0' - 9 k ) j F% \0 U . . . ,9 n ) (2.10) 

e = e ' V fc =i / 

There is also a further equation related to bound states which we do not need in the sequel. 
These equations are supplemented by the assumption of maximum analyticity (i.e. that the form 
factors are meromorphic functions which only have the singularities prescribed by the equations) 
and possible further conditions expressing properties of the particular operator whose form factors 
are sought. 

2.2 The thermal two-point function 

Let us take a field theory at finite temperature T, which can be formulated in a periodic 
Euclidean time 

t = t + R where R = l/T (2.11) 
Our aim is to determine the following correlation function 

t)O 2 (0)) R = 1 Tr(e -^) ( 2 - 12 ) 

Naively, one can proceed by inserting two complete sets of states to obtain the spectral represen- 
tation 

J- e -(R-t)E„ e -tE me ix(P n -P m ) (q) | m j ( m p 2 ( ) | n j 

<CMM)0 2 (O))* = ™ _ REn (2.13) 

2^e 

n 

where E m and E n denote the total energies, while P m and P n denote the total momenta of the 
states inserted. Using asymptotic completeness, the sets of states {|ro)} and can be chosen 

as the basis of asymptotic in (or out) states. In this case, the matrix elements appearing in the 
above formula are just the form factors of the operators 0\ and 2 . 

However, the above expression is ill-defined because according to the crossing relation (2.5) any 
term in which there is at least one particle with the same rapidities in the states |n) and \m) both 
matrix elements contains a 8 function term, and so the expression contains squares of 5 functions. 
Similar divergences occur in the partition function in the denominator. A standard combinato- 
rial consideration of disconnected terms shows that the divergent parts cancel out between the 
numerator and denominator [42]. The issue is therefore to compute the finite remainder. 

The simplest idea is to put the system in a finite spatial volume L with periodic boundary 
conditions 

x = x + L (2.14) 
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in which case the expression becomes 



R Tr L (e-^Oi(x,t)O 2 (0)) 



where Trr, denotes the trace over the finite-volume states, Hl is the Hamiltonian in volume L. 
This expression can be expanded inserting two complete sets of states 

Tr L (e- RHL O 1 (x,t)O 2 (0)) =Y,e~ REn{L) (n\Oi(x,t)\m) L (m\O 2 (0)\n) L (2.16) 

m,n 

where the matrix elements of local operators are also taken in the finite volume system. 

2.3 Form factors in finite volume 

The next ingredient we need is the description of form factors in finite volume. Previously 
this was achieved using semi-classical techniques [55, 56], but that is not suitable for our purposes 
here. We need a formalism that gives the exact quantum form factors to all orders in L~ l (i.e. up 
to corrections that decay exponentially with the volume). The relevant results were derived by us 
in [44, 45], Following our conventions in those papers, the finite volume multi-particle states can 
be denoted 

\{h,...,I n })L (2.17) 

where the I k are momentum quantum numbers. We can order the momentum quantum numbers 
in a monotonically decreasing sequence: I\ > • • • > /„, which is just a matter of convention. The 
corresponding energy levels are determined by the Bethe-Yang equations 

e imLsinh0 fe "Q _ fy) = 1 (2.18) 

l^k 

We define the two-particle phase shift 5(9) by the relation 

S(6) = -e im (2.19) 

where the — sign ensures that (due to the generic feature S(0) = — 1 and the bootstrap relation 
S(9)S(—9) = l) 1 the phase-shift can be defined as a continuous and odd function of 0. In addition 
we introduce the following notation 

m = ^ (2.20) 

for the derivative of the phase shift. Using these definitions we can write 

Q k (e 1 ,...,e n ) = mLsinh9 k + ^2S(e k -ei)=2xI k , k = l,...,n (2.21) 

where the quantum numbers I k take integer/half-integer values for odd/even numbers of particles 
respectively. Eqns. (2.21) must be solved with respect to the particle rapidities k , where the 
energy (relative to the finite volume vacuum state) can be computed as 



cosh 4 (2.22) 



fc=i 

up to corrections which decay exponentially with L. The density of re-particle states in rapidity 
space can be calculated as 

p(e 1 ,...,9 n )=detjW , J £) = dQ k (9 1 ^..,9 n ) M = 1 _. )n (223) 



1 S(0) = — 1 is valid in any known integrable model except for the free boson, which we do not consider here 
(as the issues of this paper can be solved trivially), while the other relation is a consequence of the unitarity and 
Hermitian analyticity of the S matrix. 
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The finite volume behaviour of local matrix elements can be given as [44] 

<{/(,..., r m }|O(0, 0)\{h, . . . , I n }) L = ^ fn + i*,---A+i*,Oi,---,On) + 0{e -, L) (2 24) 



8' 



where 8 k (0' k ) are the solutions of the Bethe-Yang equations (2.21) corresponding to the state 
with the specified quantum numbers I±, . . . ,I n . . . ,I' n ) at the given volume L. The above 
relation is valid provided there are no disconnected terms i.e. the left and the right states do not 
contain particles with the same rapidity, i.e. the sets . . . , n | and . . . , are disjoint. 

The lower limit on the exponent p is independent of the operator and the states considered (it is 
related to the bound state pole structure of the infinite volume scattering theory). 

It is easy to see that in the presence of nontrivial scattering there are only two cases when 
exact equality of (at least some of) the rapidities can occur [45]: 

1. The two states are identical, i.e. n = m and 

{l{,...,I , m } = {h,...,I n } (2.25) 

in which case the corresponding diagonal matrix element can be written as a sum over all 
bipartite divisions of the set of the n particles involved (including the trivial ones when A 
is the empty set or the complete set {1, . . . , n}) 

({h . . . I n }\0\{h . . . I n }) L = ^ AC{1 ' 2 '- n} (f \ m ' + 0(e-" L ) (2.26) 

p{{l,...,n}) L 

where \A\ denotes the cardinal number (number of elements) of the set A 

p({k 1 ,...,k r }) L = p(8 kl ,...J kr ) (2.27) 

is the r-particle Bethe-Yang Jacobi determinant (2.23) involving only the r-element subset 
1 < k\ < ■ ■ ■ < k r < n of the n particles, and 

■F({fcl,...,*r})i = Fi r (§ kl ,...,e kr ) 

F$ l {d u ... i d l ) = UmFglPi +iir + e,...,Q 1 +in + e,6 1 ,...,e l ) (2.28) 

is the so-called symmetric evaluation of diagonal multi-particle matrix elements. 

2. Both states are parity symmetric states in the spin zero sector, i.e. 

{Ii,...,I n } = {—!„,... 

{/{,..., 4} = {-4,...,-/(} (2.29) 

Furthermore, both states must contain one (or possibly more, in a theory with more than 
one species) particle of quantum number 0, whose rapidity is then exactly for any value of 
the volume L due to the symmetric assignment of quantum numbers. Writing m = 2k + 1 
and n = 21 + 1 and defining 

Tk,l{8\, - ■ ■ , S' k \0i, ... ,9i) = 

lim F 2k+ 2i+2(i'K + 0[ + e, . . . , m + 6' k + e, m — 6' k + e, . . . , m — 6[ + e, 
i7r + e,O,0i,...,0,,-0,,...,-0i) (2.30) 
the formula for the finite-volume matrix element takes the form 

({/[, ...,4,0, -I' k , ...,-l[}\0\{h,..., I h 0, , -h}) L (2.31) 

/ ~ - ~ \ -1/2 

= ( P2k+i(6[, ... ,9 k ,0, —6' k , • • • , -0'i)p2i+i(6i, ■■■ ,0i,0, —0i, . . . , —0i)J 
x Fk,i{0'i, ' ' ' , jfc I 1 3 * * * ? 

+mL F 2fc+ 2/(ivr + §[,... ,itt + 0~' k ,m - 0~' k , . . . , iir - 6[, 0~i,...,0 h -0 h ...,-8] 
+0(e-^ L ) 
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2.4 The form factor expansion using finite volume regularization 

Using the finite volume description introduced in subsection 2.3 we can write 



z 

N,M 



where 



C NM = E E ({h...I N }\Oi(0)\{J 1 ...J M }) L x 

Ii...In J\---Jm 

{{J\ ■ ■ ■ JM}\O 2 (0)\{h ■ ■ ■ I N }) L e i{Pl - P2)x e- E ^ R -^e- E ^ (2.33) 

and Ei 2 an d -Pi, 2 are the total energies and momenta of the multi-particle states . . . In})l 
and |{ Ji . . . Jm})- The task is to calculate the sum in finite volume and then take the limit 
L — > oo. 

It is easy to see, that the terms N = 0, M = . . . oo and N = . . . oo, M = add up to 
zero-temperature correlation functions 



M=0 N=0 



lim V Com = (Oi(x,t)O 2 (0)) lim V C N0 ) = (Oi(x,R-t)O 2 (0)) (2.34) 

A 

In particular 



lim CW = ^ / ^...^^(fc,...,**)!^^^^^ 

L->oo Ml J Z7T Z7T 

(2.35) 

and similarly for Cjvo- 

Let us introduce two auxiliary variables u and v to keep track of the orders of e~ mt and 
e - m (R-t) ^ at the end both wilJ be get tQ ^ Then (232) takes the form 

(Oi(x,t)O 2 (0)) R = | E u N v M C NM (2.36) 

N,M 

We define a similar expansion for the partition function 

Z = ^(uv) N Z N (2.37) 

N 

with Zjy denoting the iV-particle contribution to the partition function. The first few terms are 
given by 

Z = l Z l = Y,z~ ElR Z 2 = Y,z- {Ei+Ej)R (2-38) 

The inverse of the partition function is expanded as 

Z" 1 = J2(uv) N Z N (2.39) 

N 

where 

Z = l Zi = -Zi Z 2 = Z\ - Z 2 (2.40) 
Putting this together we can rewrite the expansion as 

(0i(M)O 2 (O)>£ = E uN v N£) NM (2.41) 

with 

£>nm = ^2 Cn-i,m-i z i (2.42) 
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The first few nontrivial terms are given by 

Dim = C\m - ZxCq^m-i 



(2.43) 



D 2 M = C2M - Z\C\ t M-l + {Z\ — Z 2 )Cq } M-2 

In this way we produce a double series expansions in powers of the variables e~ mt and e - m ( fi_t ) . 
Since these variables are independent, each quantity D^m must have a well-defined L — > 00 limit 
which we denote as 

D NM = lim D NM (2.44) 

L— »oo 

and we obtain that 

<0i(x,t)0 2 (O)> fl = lim (d(x,t)C» 2 (0))f = V D NM (2.45) 

L— >0O * — ' 

N,M 

The reordering of the series in (2.41) using the coefficients Dnm is also an integral part of Essler 
and Konik's calculation in [53]; we used a similar reordering for the expansion of the one-point 
function in powers of e~ mR [45]. 

Note that individual terms contributing in (2.42) to Dnm may contain divergent pieces which 
scale with positive powers of L. Similarly to the considerations for the one-point function in [45], 
it turns out that the TV-particle terms which are most singular in the large- L limit carry a factor of 
(mLe~ mR ) N . Since it is also necessary that the exponential corrections to the finite volume form 
factors in eqns. (2.24) and (2.26) are small, the finite-volume expansion is valid in the domain 

l«mi« e mR (2.46) 

For the limit (2.44) to exist the positive powers of L must drop out, therefore one can understand 
the L — > 00 limit of the series as an analytic continuation to very large values of L outside the 
domain (2.46). Eventually, the condition that the coefficients D nm must have a finite large volume 
limit can be used as a nontrivial check to verify our calculations. 

It is evident from (2.33) that the quantities Dnm with N > M can be obtained from those 
with N < M after a trivial exchange of t and R — t. Therefore we will only consider the case 
N < M. 



3 Warm-up example: the zero-temperature three-point function 

Before tackling the central issue of the paper, we consider a simpler problem which allows us 
to introduce the central ideas without too many complications. Let us consider the three-point 
function 

(0\O 1 (t l ,x l )O 2 {t 2 ,x 2 )O 3 {0)\0) (3.1) 

in the Euclidean theory with non-compact time direction (i.e. T = 0). Suppose that t\ > t 2 and 
to shorten the formulae we also omit the dependence on x\ and x 2 (it can be reintroduced easily). 
The spectral decomposition takes the form 

(0|Oi(t 1 ,n)O 2 («2,a;2)O(0,0)|0) 

= ^(0|Oi(0)|m)(m|O 2 (0)|n)(n|O 3 (0)|0)e- Sm( ' 1 -* 2) e- E " t2 (3.2) 

m,n 

What makes this example simpler is that the disconnected terms appear linearly since they only 
enter from the 2 matrix element. Following the example of the three point function, we can 
introduce a finite volume regularization 

(0\O 1 (t 1 ,x 1 )O 2 {t 2 ,x 2 )O(0,0)\0) L = (3.3) 

N,M 
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where 

,(3) 



C nm = E E (0\O 1 (0)\{I 1 ...I N }) L ({h...I N }\O 2 (0)\{J 1 ...J M }) L (3.4) 
Ii...In Ji-.-Jm 

x ({Jx... J M }|O 3 (0)|) L e- £l(i)(il - i2) e- S2(L)t2 (3.5) 

There is no denominator Z to supply counter terms for the L dependence, therefore each of these 
expressions must have a finite limit as L — > oo: 

= lim Cfl, (3.6) 

L— >oo 

and 

<0|Oi(ti,xi)O2(* 2 ,x 2 )O(0,0)|0) = E ( 3 - 7 ) 
Eventually, it is trivial to write down some terms of the expansion: 

vf M = (o 1 )^J^...J^F^(e u ...,e M )F^(OM + i^...,e 1 + ^) 
( M \ 

x exp —mt2 cosh 9j 

( N \ 

x exp —m(ti — t2) J ^cosh0 i (C 3 ) 



ITT) 



(3.« 



and also 



£>§> = f ™1 f ^ F ^ F ^ ( 1 + ^^ )F 03 e -rn( t2 - t Ocosh eie - mil cosh i 

since the two-particle form factor has no kinematical singularities. 



(3.9) 



3.1 The contribution 

The first nontrivial contribution is , for which the finite volume expression is 
n (3) _ sr^ F 1 01 F 3 c ' 2 (6'i + in, 9' 1 ,9' 2 )F 2 ° 3 (9' 2 , 9[) - m (ti-t a ) cosh^ „-mt 2 (coshe;+coshe;) /o-,m 

12 ~ n (ft\ n (ft' ft'\ (J.1UJ 

where the rapidities satisfy the appropriate Bethe-Yang quantization relations. We can make a 
choice whether to perform first the one-particle or two-particle summation. The latter proceeds 
by an application of the multi-dimensional residue theorem (A. 2) and illustrates one of the central 
ideas that make the expansion of the thermal correlator feasible. 

3.1.1 Summing over one-particle states first 

We can substitute the sum over I\ 

V^W ddl Pl{ ° l) (3111 
h h ■ 7C/ i 

where 

Qi(0i) = mLsinh^ p^) = Q'(9 1 ) (3.12) 
and Cj 1 are small closed curves surrounding the solution of 

Qi(0i) = 2ttJi (3.13) 
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in the complex 9\ plane. Now we open these circles and join them to obtain the contour 

C = C+ + C_ (3.14) 

the C+ running from oo + ie to —oo + ie (i.e. backward in 9\) while C_ running from —oo — ie 
to +00 — ie (forward in Ke#i). However, by this operation we also include the contribution of 
two poles at 0\ = 9[ and 9\ = 9' 2 which must be subtracted. Using (2.10), the singularity of the 
integrand at 9\ ~ 9[ can be written as 

1 iF^{l - S(6[ — 9' 2 ))F^ >2 F^ )a (6' 2 ,9' 1 ) mtl CO sh6>i -mt 2 cosh 6^ /n i z\ 

p 2 (6' 1 ,9' 2 )(e i Q^)-l) 1 ' 

The quantization relation of the two-particle state can also be written as 

e imLsmh9' lS ^_ e ^ = 1 (3 ^ 

and so we can rewrite (3.15) as 

1 77^2 77C3 1 Ql a' \ 

1 lr l r l c 2 \ a li u 2) -mil cosh 6[ -mt 2 cosh 6' 2 / n 1 7 \ 

9 1 -e[ P2 (9[,9' 2 ) [6 - U) 

The contribution of this singularity can be evaluated as 

d6i 1 iF^ )l F^ 2 F 2 s (9' 1 ,9 2 ) mt coshe / - m t 2 cosh 6^ 
e i e * — 



,271-01-01 P2(0' 1 ,9' 2 ) 

rpOi rp0 2 rpOs (nl nl \ 

r l r l r 2 \ u l' u 2) -mix cosh 9[ -mt 2 cosh 6' 2 

P2(9[,9' 2 ) 



(3.18) 



The contribution of the 9\ = 9' 2 pole can be obtained in a similar way. These must be subtracted 
from the 9\ integral and therefore we obtain 



d9i F^ )l F^ 2 (9i + iTT,9' 1 ,9' 2 )F 2 ' i (9 2 ,9' 1 ) m{ti cosh di mt2 (cosh e , +cosh ; 
c 2vr p 2 (9[,e' 2 )(e^)-l) 



Ji,J 2 

771O1 7-1O2 rpO:i I nl nl \ 
I r l r l r 2 \ a l' u 2> -mti cosh 0'-mt 2 cosh 6^ 

P2{9' l A) 

rpOi T?0 2 77O3 (nl nl \ 
I r l r l r 2 \°2-> u l) -mt-i cosh 6' 2 -mt 2 cosh 6^ 



(3.19) 



Taking the large L limit, we can substitute the discrete sum with an integral 

Jl<J2 J J 



and so we 


obtain 






r (3) 

°12 


1 rrd9' 1 d9' 2 

~2jJ ^27 








00 dBx F^F^ 2 (9 


1 +i(vr + e),( 


^,0 2 )if 3 (0 2 ,0i) 




-oo 27r 


gi<2l(0l+ie) 


- 1 




00 d0! F° 1 F° 2 (9 


1 +i(7r - e),( 


?i,^)lf 3 (Ml) 


-00 2ti" 


e iQl(0i— ie) 


- 1 



_|_ jt'Oi JT 1 ^ p®3 f^Ql Ql ^g— mil cosh 0^ —mt 2 cosh 0j 
_|_^7iOi j?0 2 pOz ^gl qI ^g— mil coshSj — m*2 cosh 0j 



g— m(ti — t 2 ) cosh(0i+ie)— mt2 (cosh 6^+cosh 6*j) 



-m(ii — t 2 ) cosh(0i —ie)—mt 2 (cosh 0^+cosh 8' 2 ) 



(3.21) 
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Note that 

iQl (#1 i *e) = ^fmL cosh #1 sin e + imL sinh 0\ cos e 
and therefore the L — > oo limit yields 



(3.22) 



(3) 
12 



„ — m(ti —42) cosh(9i+ie)— m42(cosh 6^+cosh $2) 
j^pO\ p02 pOi ^qI qI ^g— mti cosh 0' 1 —mt2 cosh0' 2 

■j-F® 1 F® 2 F® 3 (#2 ^^e -7 ™* 1 coshe 2 — m *2 cosh^ 



(3.23) 



3.1.2 Evaluating the two-particle summation first 

Using the multi-dimensional residue theorem (A. 2) we can represent the two-particle sum as 

1 ^ f f d6[ d0' 2 1 

f CjiJ2 2vr 2vr ( e iQi(o[,e 2 ) + i) ( e iQ2(M 2 ) + i) 



E 

Jl>J2 



E 

Jl>J2 



(3.24) 



where Cj x j 2 is a multi-contour (a direct product of two curves in the variables 6[ and 2 ) sur- 
rounding the solution of 



where due to the definition 



mLsinh^i + J(0£ - 6» 2 ) = 2vrJi 
mL sinh 6» 2 + <5(0 2 - 8[) = 2vrJ 2 

5 = -e iS 



(3.25) 



Ji and J2 take half-integer values. Since the form factors vanish when any two of their arguments 
coincide, we can extend the sum by adding the diagonal 



E 

Ji>Ji 



IE 



and so 



where 



C 



(3) 



12 



E 

h 



Ci 2 (gQ 



(3.26) 



E 



Ji,J 2 



O 2tt 2tt (e'ft( 9 ;.«2) + 1) (^2(^2) + 1) 



xe 



-m(ti — tz) cosh 6>i — mt2 (cosh 0j+cosh #2) 



(3.27) 



Now we open the contours to surround the whole of the real 0[ and 0' 2 axes (but close enough so 
as to avoid all singularities of the S matrix). Just as before it is necessary to subtract the contri- 
butions of any singularities encountered in the process. There are two classes of such singularities: 



• 0[ = 61 and e iQi + 1 = 

• 0*2 = 01 and e iQl + 1 = 

There are no triple singularities because 0\ satisfies 



?7iLsinh#i = 2irli 



h G 



(3.28) 
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and it is impossible for the three Bethe-Yang conditions Qi, Q 2 and Q3 to be satisfied simultane- 
ously. As before, it is enough to evaluate the first case; the second can be obtained by swapping 
9[ and 9' 2 . In the large L limit 



JQl,2(e' 1 ±ie 1 ,e' 2 ±ie 2 ) 




(3.29) 



therefore we obtain 
612 = 




d0' x d9' 2 



C++ 



2vr 2tt 



Ff 1 F° 2 (8i + ITT, 9[, 9 2 )F 2 ° 3 (9' 2 , ^Je-TOCti-taJooahfli-mtaCooBh^J+coA^a)) 



Res 



J 2 GZ+l/2 Q2=J2 



( e iQi(0i,0S) + 1) (^M) + !) 



xe 



-m(ti — ta) cosh 6<i —mt2 (cosh 9^ +cosh 9' 2 ) 



(3.30) 



where C++ denotes the part of the two-particle multi-contour that survives in the L->oo limit; 
it corresponds to an integration parallel to the real axes in 9[ and 9' 2 with a shift in the positive 
imaginary direction, i.e. 



Recalling (3.25) we get 



Res 



F° 2 {Ox + ivr, 9[ , 6> 2 )F 2 ° 3 (9' 2 , 0J ) e - m (* 1 cosh ^ " m * 2 ( cosh °i +«»h^) 
pOTCepD + 1) (e^W^) + 1) 



-mti cosh 61 —mt2 cosh 



(1 - S(0i - 2 )) (-i) (mLcosh0 2 + <p(9' 2 - 9{)) 

(3) 

Substituting these into the expression for C l2 and taking L->oowe obtain 



(3.32) 



D 



(3) 
12 



/CO Jfl 
f^if 3 if 2 (<? 1 + *7r,ei,^) J Ff H^^i) (3-33) 



-m(ti — £2) cosh 9\ —mti (cosh(0^ )+cosh(0 2 )) 



+ F^F^F 2 °\9[,9 / 2 )e 
+ F^ 1 F^ 2 F 2 ° a (9' 2 ,9' 1 



mt2 cosh (9j g— mti cosh 2 
g— m42 cosh 6>2 g — m *l cosh 6^ 



After substituting 0^ 2 — ^ ~~ 2 an< ^ using 

F^(-9' 1 ,-9' 2 ) = F^(9 2 ,'9' 1 ) 
we can make a combined shift of the three contours in the triple integral term to obtain 



D 



(3) 
12 



d0j d0 2 
2vr 2vr 



^ pOi P C 2 / fl 



+ z( 7 r- e ),0i ) 2 )F 2 O3 (0 2 ,0i; 



(3.34) 



^ g— m(ti — ta) cosh(6i — mti (cosh 0^+cosh 8' 2 ) 

+ J?® 2 (J) 1 ff f ) g — m * 2 cos ^ g — mi 1 cos ' 1 ^2 



+ F^F^F 2 ° 3 (9' 2 ,9' 1 )e 



-mt2 cosh 2 p — m t 1 cosh 6^ 
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It can easily be shown that this expression agrees with (3.23); the difference due to the e 
change drops out: 

-m(t2 — ti) cosh(0i+ie) p — mti (cosh 9' 2 ) 



-m(t2 —ti) cosh(Si —ie)^—mt\ (cosh 6^+cosh 6' 2 ) 



/oo in 

/oo in 
^*f 1 *f a Vi + K* - e),9[,9 2 )F 2 D ^9' 2 ,9' 1 )e 

= i ( Res + Res J Ff^Ff 2 ^ + in, 9[, 9' 2 )F° 3 (9' 2 , ^)e^^-*0«^^ e - mtl t c08h ^ +OOBh ^ ) 
= - (ifVf 2 (Vf 1 (0 2 , (9i) - F 2 Gl (0j, 0£)) e" m * 2 cosh0 'i e - mi i coshe 2 _ ^,)) (3.35) 
and the integral of the last expression over 9[ 2 vanishes by symmetry 

3.2 D® 

The new aspect in this case is that the contribution must be split into two parts: one with the 
two two-particle states being different and the "diagonal" when these states are the same. Using 
the finite volume form factor formulae (2.24) and (2.26) we can write 



(3) _ ^ if 1 (e l ,e 2 )Ff i (9 2 + iTT, 01 + in, g 3 , g 4 )if 3 (9' 2 , 9[] 

P2(e 1 ,e 2 ) P2 (e' 1 ,9' 2 ) 



^ g— m{t\ — t2)(cosh 6>i+cosh 62 )g— m<2(cosh 6^+cosh 0' 2 ) 

f 4 °=(«i.»2) + Eg* ( w («o + Pl (s 2 )) + l&mlthM 



+ £ 



P 2 (01,0 2 ) 2 

x F 2 Vl (9 1 ,9 2 )F 2 :)3 {9 2 ,9 1 )e- mtl(coshei+coshe2) (3.36) 
The diagonal part can be written as 



^3)dia g = yy 
{h,h} 



F^(9 1 ,9 2 ) + F 2 °J (pi(fli) +pi(fl 2 )) + (0 2 )p 2 (0i,0 2 ) 
p 2 (0i,0 2 ) 2 



x F 2 Vl (9 1 ,9 2 )F 2 ° 3 {9 2 ,9 1 )e- mtl{coshei+coshd2) 

t^oo \ //^^^H^i^ 2 )(0 2 )^ 3 (^ 2 ^i)e- m ' l(coshei+coshe2) (3.37) 
For the non-diagonal part we need to evaluate 

= 1 E Ff^ + i, e^ M ) K tutMM) (3.38) 
2 f-f P2{9i,9 2 ) 

where 

i^T tl ,i 2 (<9i , ^ 2 , , ^2 ) = ^2° 3 (^1 ? ^2 )i 71 2 ° 1 (^i ^ ^2 )e~ m( * 1 (cosh 1+ cosh 2 ) e -^t 2 (cosh 0^ +cosh ^ ) (3,39) 

(again we extended the I\ < I 2 summation by symmetry and included the diagonal I\ = I 2 where 
the form factors vanish). Using the residue trick it can be represented as 

1 ^ r l d&\ d9 2 Ff 2 (0 2 + in, 9 X + in, 9[ , 9' 2 



where 



2^ft T 2^ ( e «W«iA> + I) (e *Q 2 (^2) + i) Kwih'WM (3.40) 



Qi(0i,0 2 ) = mLsmh9 l + S(9 1 - 9 2 ) 

Q 2 (9i,9 2 ) = mLsmh9 2 + 6(9 2 - 9 1 ) (3.41) 

To open the contour we need to find the singularities that do not result as solutions of Qi t 2 = 
2nl\ 2. There are the following possibilities: 
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• Q 2 = 2ttI 2 and 6>i = 9[ or Q x = 9' 2 

• Qi = 2vr/i and 9 2 = 9[ or 9 2 = 9' 2 

• 9\ = 9[ and 9 2 = 9' 2 or 9 2 = 9[ and Q\ = 9 2 . Albeit the form factor is regular at this point, 
the denominator has a double zero due to the quantization condition satisfied by 9[ and 9' 2 : 

mL sinh 9[ + 5{9[ - 9' 2 ) = 2vr J x 

mLamh9' 2 + S(9 2 -9' 1 ) = 2vrJ 2 (3.42) 

These were excluded and their contribution calculated in the diagonal part C 22 ^ g . 

As an example we consider the contribution from the singularity Q 2 = 2ttI 2 and 9\ = 9[. Using 
(3.42) we can evaluate 

e iQi(fli,fc>) = _ e imL sinh 6[ s ^ _ = _ S r^/ _ e' l )S{9' l - 9 2 ) (3.43) 
and the appropriate residue takes the form 



1 , 2 i(l - S(9 2 - 9' 1 )S(9' 1 - e f 2 ))Ff*(0 2 + iTr, 9' 2 ) 

2 (1 - S(9' 2 - 9' 1 )S(9[ - 9 2 ))(-i)(mL cosh 9 2 + tp(9 2 - q'^ 1 ^ 1 ^ 2 ^ 1 ^ 

When we sum over I 2 we must exclude I 2 = J 2 i.e. the term 9\ = 9[ and 9 2 = 9' 2 (this singularity 
was taken into account in C^ dia " s ). However, its contribution to the I 2 sum is 

1 F°' 2 

T^K^tMAAX) (3-45) 



2 (mLcosh# 2 + (p(9' 2 - 9[)) 

which vanishes when L — > oo. 

One can calculate the contribution of all other singularities in a similar way. Taking L — > oo 
the final result is 

*>% = III ~ II P ^§Ff 2 ^ + -,^ 1+ ,vr,^,^(^ 1 ,^3 ( ^,^ 



+ 



2 J Jc++ 2vr 2vr J J 2vr 2vr 

^ g — m(ti— 12 ) (cosh 0i +cosh 02 ) g — m *2 (cosh 0^ +cosh 2 ) 
^2 

2vr 

X g — m(ti— 12 ) cosh (?2 g — mt 1 cosh 0^ ^— mt2 cosh 2 

1 /" f d9id9 2 ^Q xfQ Q ^, n „\ w 3( n n s -mti (cosh fli+coshfla) 



^^i^i(e 1 ,e 2 )(<D 2 )F 2 U3 (e 2 ,ei)e^ 1 ^^ +cosa ^ (3.46) 

2 J J 27T Z7T 

where the (9\,9 2 ) contour C++ is specified in (3.31) 



4 Evaluating the thermal correlator 

Now we show how evaluate the series (2.45). Since according to (2.34) the contributions Don 
and Dmq are identical to terms contributing to the zero-temperature two-point function, the first 
nontrivial temperature correction is given by Dn, which is evaluated in the next subsection. The 
contributions -D12, D\ n for arbitrary n > 2, and D 22 are calculated in subsections 4.2, 4.3 and 4.4 
respectively. The final expressions (which constitute the main results of this work) are given by 
equations (4.8), (4.22), (4.47) and (4.89), respectively. 
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4.1 The Du correction 

According to (2.43) and (2.44) 



Du = lim Du 

D u = C n - ZxCoo (4.1) 

where 

Coo = (Ox) (0 2 ) (4.2) 

and 

C u = £ ({/}|O 1 (0)|{J}) L {{J}\O 2 (0)\{I})l e^-^e-^-'h-^ (4.3) 
i,j 

where E\, E 2 and p±, p 2 are the finite size energies and momenta of the one-particle states. Using 
the Bethe-Yang quantization conditions (2.21) we have 

rnL sinh 9 = 2ttI , mL sinh 9' = 2ir J (4.4) 

and 

E\ = m cosh 9 , p\ = m sinh 9 
E 2 = m cosh 9' , p 2 = m sinh 9' 

According to (2.24) and (2.26), the two-particle matrix elements are given by 

<{I}|Oi(0)|{J}> £ = ^== + ^W (4-5) 

({J}|O 2 (0)|{/}) i = ^ { *t iir f £ +8ij(0 2 ) (4.6) 



Substituting the above formulas into (4.3) one obtains 

Z7T,fl')F 2 02 (i 

W)pW) 



I, J 

^ 1 (9 + i n,9) _ EiR 

pi (00 t^/z^ Pl{6) 



+ ^(0 1 )(0 2 )e- £lfl (4.7) 



The last term in (4.7) is O(L), but it is canceled in Dil by the term Z\Cqq. All the other terms 
have a finite limit as L — > oo which can be written in the form 

Dn = J M. j ^ F o^ e + i ^^^ F o 2 ^, + i ^^y(^e-M) m ^ 

+ ((d) F$> + {0 2 } Fg) J ^e-mRcos^e (4 8) 

Note that according to (2.10) the two-particle form factor does not have kinematical singularities. 
The result (4.8) was first obtained in [57]. 
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4.2 More than just a warm-up: D 12 

According to (2.43) and (2.44) 

D12 = lim (C12 - ZiCoi) 

where 

C ^ = EE <{i}|Oi(0)|{Ji, J 2 }> L x ({J l5 J 2 }|O 2 (0)|{/}) L e^- p2 ^e- El ^)e-^* 
/ J1J2 

V- F 3 ° 1 (g 1 +»7r,gl,^)if 2 (g 1 +z7r,^,gl) , (4 ' 9) 

" V ^ ( 11 1! 2) 

J£ (Q^ Ql Ql \ gimx(sinh 9i — sinh 6[ — sinhSj)^— m(R— t) cosh (?i ^— mt(cosh 0j+cosh Sj) (4 10) 

7 _ / \^ a -ER \ ( f dQ J? 1 rp0 2a -imxshe-rntche\ 

ZiCi-l^e I (J -F, e j 



and 



4.2.1 First summation: one-particle states 

We first perform the summation over /. The quantization condition reads 

Qi(6i) = mLsinhOi = 2irl p\ = = mL cosh 9\ 

do i 

with IeN. Therefore it is possible to convert the summation into a sum over contour integrals 
d0 1 F^(9 1 +m,9' 1 ,9' 2 )F^(9 1 +m,9' 2 ,9' 1 ) , R) 



fi^f 2^ ^e[A) (MM&M&TTi (4 ' n) 

Jl J2 i 

In order to open up the contours one has to calculate the surplus singularities of the integrand, 
which appear at 9\ — > 9[ and at 9\ — > # 2 - Each of the form factors have first order poles, therefore 
the singularity is a second order pole. In the following we calculate the residue at 9\ — > 9[; the 
case 9\ — > 9 2 will be given by a change of variables. 
The residue of the form factors for 9\ — > 9[ read 

if 1 (9 1 + in, 9[,9' 2 ) = i(l - 5((9i - # 2 )) + . . . (4.12) 

Let us introduce the connected part of the three-particle form factor as 

^MWiA) = g Vm (if ^ +in,e' 1 ,e' 2 ) -i(i -s(e[ -^)) J£l) (4.13) 

The connected form factor defined above still has a pole at 9[ = 9' 2 . In fact, the singularity 
structure of the original form factor near 9\ = 9[ = 9' 2 is given by 

if ^ = 2iif 1 (^-Ljr ~ o^qt) (4-14) 

and after subtracting the first pole there remains the second one leading to 

F^MK^) = -2iF° 1 ¥ ± ¥ + ... (4.15) 

Also, it can be proven that 

fZWiWiA) = s(9[ - e' 2 )Fii(e' 2 \e 2 ,9[) (4.16) 
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In the case of the crossed form factor one has 



F o 2 



F^(9 l+ in,9 2 ,9[) = - S(9' 2 - ^-L_ + S(9' 2 - e^MKA) + ... (4.17) 
With these notations the residue of (4.11) at 9\ = 9[ is expressed as 
K t R J ( 6 i ^1^2) { (i S ( 9 'i ~ °2) - l)(»m ch 9[x - m sh 9[ (R-t))+ imL ch 9'^j F® 1 if 2 

+ iF^MK^F? 2 + iFgiio'iK e' 2 )F^ } 



(4.18) 



There is a similar contribution at 9\ = 9' 2 , with the role of 9[ and 9' 2 exchanged. After integrating 
over 9[,9' 2 one could make a change of variables to obtain the same contribution twice. However, 
one has to keep both residues separately because of the poles of the quantities F^ sc . Making the 
change of variables only in the regular terms one obtains the two contributions 



dsing FF =2kQ ) (9' 1 ,9' 1 ,9' 2 )((S(9' 1 -9' 2 ) - l)(imch9[x - msh9' 1 (R - t)) + imL ch 9[) F® 1 F° 2 
ssing FF = l K^(9' 1 AA) (F^MKe'jF?* +Fgl(9' l \9' l ,9' 2 )F? 1 ) + 

iKl R J(e' 2 ,e' 2 A) (f^MKo'i)f^ + Fgi(e' 2 \e' 2 ,9' 1 )F^) 

(4.19) 



Now is is possible to perform the summations over 9[,9 2 . The O(L) term of (4.19) can be 
transformed in the L — > 00 limit into 

mLF 1F 2 Q gg-^A^-mchfljA Q g. ch0 / e -mch^ _ 

/jnl (4.20) 
—im sh 8i,x 
2n 

The first term gets exactly canceled by Z\Cq\ leaving only the finite contribution 

— F 0l F° 2 t < ^k e - mche 2 ( , R + t )- irnshe 2 x (4 21) 

1 1 J 2-7T 

The remaining terms of (4.19) are regular, therefore it is allowed to replace the summation over 
9[,9' 2 with the appropriate integral. The final result is 

+ / 5||i^?(0iA^)m 

+ \ J ff^^i^^^^^l^,^ + Ft(0[\e[,9' 2 )F^) + (9[ o 9' 2 )} 



- F® 1 F® 2 I d ® 2 c ~mch 9!,(R+t) —im shdLx 



(4.22) 



4.2.2 Performing the two-particle summation first 

We can express C12 as 

sr^_srll d0 i d °2 F3 1 (Oi + «r, 6[ , 9> 2 )F° 2 (9' 2 + in, 9[ + in, 9,) {R) 
1 Pl{9l) fa, J ^2 2n 2n (e^'M + 1) ( e «M*i^) + 1) l ' x { h 11 2) 

(4.23) 



17 



jy~(R)(a Ql a 1 \ imx(smhdi— sinhS^ — sinh8' 2 ) — m{R— t) cosh 61 — mt(cosh (9j+cosh Sj) 

Now we open the multi-contour to surround the real axes in 8[ and 9' 2 ; however, we encounter 
some "surplus" singularities: 

• QF poles, where the singularity in one of the variables arise from a Q-denominator while in 
the other from a form factor: 

e[ = e 1 , e^'^) + i = 
e' 2 = e x , e lQ i'(<^ + i = o 

• FF poles, where the singularity in both variables comes from the form factors: 

8 l = 81 , 6*2 = 6*i 

(note that positions where there is only a singularity in one of the variables do not contribute, 
as the contour in the other variable can be shrunk to a point). In the following we calculate the 
contributions of these singularities. 

The FF singularity 

We can write 

if 1 (8, + zvr, 9[,8' 2 ) = —^—t (1 - S(8[ - 8' 2 )) if* + {S(8' 2 - 9[) - l) if 1 

y i - y i u i- 8 2 (4.24) 

where F 3c l is the regular part of the form factor around the singularity. The pole contribution is 
then 

1 r d8[ r d8> 

2 J Ql 2vr f 9l 2tt (e^i'M) + 1) (jQAW + i) ^ ^ 



X 



(1 - - 8' 2 )) if 1 + ^l^i (fif(6K 2 - ei) - 1) Ff 1 + f£(fliKy 2 )) 
(1 - S(8[ - 8' 2 )) if 2 + ^-L-i (S(8' 2 - 8[) - 1) if" + 6K 2 )) 



i^cc does not contribute since then either the 9[ or the # 2 integration contour can be contracted 
to a point. For similar reasons, the only terms that could give a nonzero contribution are the 
"cross-terms" 

1 / dB[ f d8' 2 , , 

2 f 9l 2tt f 9l 2n (e^i'M) + 1) (jQAo'itt) + l) { 1 2j 

x {h=of (1 " 5(ei " 6 ' 2)) F ^ x ^ " 9[) " 1} Fl ° 2 ) 

We need the residue at 8' x = 8' 2 = 8\. Then 

S(8[ - 8' 2 ) = 5(0) = -1 

and 

e iQ 1 / (6*1,61) _ e iQ 2 /(9i,9i) _ e imL sinhfli _ ^ 

therefore the contribution of the pole is given by 

T-iOi 771O2 „— imx s'mh9i— m(R+t) cosh^i 
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After performing the 9\ sum converted to an integral we find 

Cjpp _ J ^ pOi p@2 ^—imx sinhc9i— m.(R+i) cosh 0\ ^ 2'Zi) 

which correctly reproduces the last term of D\ 2 (4.22). 
The QF pole at 0[ = 9 X 
The singular contribution is 

i r d9[ i de> 2 (gi + in, 9' 1 ,9' 2 )f^(9 1 + ^ e' 2 ,9[) rAR) , 



2 J 0l 2tt J Cj2l 2n ( e «3i'(»i.^) + l) (e^'iW + i) K ^ ^'^i'^) ( 4 - 2 6) 
where Cj 2 / surrounds the 6' 2 solution of 

Q 2 ,{9 U 9' 2 ) =mLsmke' 2 + 8{e' 2 -e x ) = 2ttJ 2 (4.27) 

where 9\ is the solution of 

Qi(0i) = mL sinhtfi = 2tt/ (4.28) 

The behaviours of the form factors are given by (4.13) and (4.17). We can separate the integrand 
into two terms according to the order of the 9[ = 0\ singularity. The first order term has the form 

9 i __i I d6' 2 i (tJ(gr _ gl , 1 , 

QF ~ 5 f 01 2vr J C]2l 2vr (e^M) + 1) (e^'M) + l) - X 11 2 lj > 

x (if 2 j# | e[ ,9' 2 ) + f^ *£• (9^,9, ; 

This can be easily evaluated: 



where 

PaWi) = -^rQAQiA) = mLcaah&z + V (9' 2 - 9 X ) 



is the density of # 2 solutions for a given #i. This gives 

S QF = -\H ^^e^ 8 ^ ~ 6l) ft-M) + F^F^(9 1 \9 1 ,9' 2 )) (4.29) 

Together with a similar term Sq F obtained by swapping 9[ 02 this gives the contribution 

K£\e' 1 ,9' 1 ,9' 2 )S(9' 2 -9[) (F^F°^9[K0' 2 ) + F°>F£(9[\9[,e' 2 ) 





d9[ r 


d9' 2 y 


II 


2vr J 


2vr . 




^9' 2 ) 





(4.30) 



The second order term reads 

*(«)/ 



D 



l 

QF 



d 9[ r d9' 2 K^ x , (9 1 ,e[,9' 2 ) 



2 J 01 2tt J Cj2l 2tt (fSQi'VM + 1) ( e ^2'M) + x ) 

x (l-^-4))(^-^)-l) F Q Q 2 
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The contribution of the double pole can be evaluated by taking the derivative with respect to 9[. 
The result reads 



D'qf 



1 



d9 2 F 1 ^F 1 ° 2 KQ ) (9 1 ,9 1 ,9' 2 ) 



2 Ic J21 2vr 



9' 2 ) + 1) 

+imLcosh^i - vp(9i - 9' 2 )S(9' 2 - 0{) 



(—imx cosh 9\ — ?7iisinh#i) (S(9' 2 — 9±) — l) 



1 

+ 2 



d9'F^F^KQ\9 1 ,9 1 ,9' 2 ) 



Cj 2 i 



2vr 



{S{9' 2 - 9±) - 1) ip{9' 2 - 9,) 



The last term with the double pole yields zero for L — > oo since it is proportional to L 2 and 
becomes L _1 after including the density. Explicitly it evaluates to 



1 f d9' F^F 1 ° 2 K^ ) (9 1 ,9 1 ,9' 2 ) 



2j Cj2l 2vr -Q' 2 ,(0i,0 3 *) 2 (02-03*) : 
where #3* is the location of the solution of 



(S(9' 2 - 0i) - 1) <p(9' 2 - 9 1 ){-S{9' 2 - 9,)) 



Q2'(e u 9 2 ) = 2irJ 2 



(4.32) 



but 



<22'(#i,#3*) 2 = (mLcosh0 3 * + <K03* -^i)) 2 = 0{L~ 2 ) (4.33) 
Even after multiplying this by the density Q' 2 , [9\ , #3* ) when converting the summation over #3* 
to integral a suppression 0{L~ l ) remains, resulting in zero large volume limit. 

Putting in the 9\ summation and a factor 2 to account for the contribution obtained by 
exchanging 9[ with 9' 2 , plus a minus sign since this is to be subtracted in the end, and adding the 
— Z\Cq\ term gives 



d ^j d ^F^{9 x M^, 



-mL cosh #1 + i P (9 1 - 9' 2 )S{9' 2 - 9{) 



(ma; cosh #1 — imtsinh#i) (S(9' 2 — #1) — l) 



m L I ^ cosh fl e - mRcosh 6 FfVf 2 / ^L e - ima:sinhe '- m * cosh 



2vr 



d9' 



2vr 



The O(L) term cancels as expected, and after a partial integration one obtains 

D Q f = - I— I ^F 1 0l F 1 ° 2 K { t f(9 1 ,9 1 ,9' 2 ){mx C osh9 1 +im(R-t)sml\9 l )(S(9' 2 
J 27T J 2tt 



End result 



Putting together everything 



2 

+ 



C++ 

1 r d9[ f d9' 2 

2 J ~2tt J ~2tt L 

(9[ o 9' 2 ) 



2vr 2vr 



+ i7r,^,ei)4?(«i^i,^) 



?0 2jr O ua i \ Q I Ql 



/ ^ F i 1F i 2K S ) ( e ^ e ^ d 2)( mxcosh0 i +im{R - t)wHx6i) {S(9' 2 - 9 X ) - l) 



2tt 
d9 



f ^_ e -m(R+t) cosh 9-imx smb. pOi p0 2 
J 27T 11 



(4.34) 
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It is a straightforward, although somewhat tedious exercise to show that the above expression can 
be transformed in the form (4.22). First of all observe, that shifting all three variables the first 
term of (4.34) can be written as 



27 



d9[ 
~2n~ 



2vr 



F^{8 1 +in,8' 1 ,8' 2 )F^ 



+ in,9' 2 ,9' 1 )K[ R) ' 



h, 6[,9 2 ) 



This differs from the corresponding term in (4.22) in the contour for 9\, which in this case runs 
below the real axis. Shifting this contour to run above the real axis one picks up the poles of the 
integrand, which can be evaluated using standard techniques. It can be shown that the resulting 
contributions are 

_ f d9[ J ^^_2 + s(e[-e' 2 ) + S(e' 2 -e[)){mx C oshe[+im(R-t)smh9' l )K { t R) (9' 1 ,9' 1 ,8' 2 ) 
J 2tt J 2ir 

(S(9' 2 -9' 1 ) - l) (f^F^{9[\9[,9' 2 ) + F^F^(8[\8[,8' 2 )) (9' 1 ,9[,9' 2 ) 



+ 



+ 



d9[ 
~2~n~ 



d8' 2 
27 



t o 0'' 



Adding these terms to the second and third lines of (4.34) one recovers (4.22). 

We wish to note, that if instead of (4.23) we had started with a similar formula including the 
factors 

1 



( e -*Qi'(0i,^) + 1) ( e -iQ 2 '(0'i,02) + l) 
we would have arrived immediately at the result (4.22). However, the calculation presented above 



is a non-trivial cross-check of our methods. 



4.3 The contribution L> ln for n > 2 



Based on the previous subsection it is now straightforward to evaluate the contribution F)\ n 
for arbitrary n. It is given by 



where 



D\n — hm (Cin — ZiCo n -i) 
L->oo 



Cl « = ^EE ({I}\O 1 (0)\{J 1 ,...,J n }) L {{J 1 ,...,J n }\O 2 (0)\{I}) L 

I J\...J n 

Xe i(Pi-P2)x e -Ei(R-t) e -E 2 t 

1 V- V- FnU d + ™A,---,0n)F n D * 1 {9 + iit,9i n ,...,9' l ) 



n\ Pl (9)p n (9' v ...,9 n ) 

^ ^J,mx(smh 8—mx y"V sinh 8j)—m(R—t) cosh 8—mt yv cosh 6j ^ 

There are additional disconnected terms in the case of n being odd, according to the rule explained 
in subsection 2.3. This happens in the presence of zero-momentum particles, which requires 1 = 
and the set {Ji, . . . , J n } to be parity symmetric; the disconnected term is given by formula (2.31). 
It is easy to show using the constrained density of states, that all contributions associated to these 
disconnected terms scales with negative powers of L, therefore we neglect them in the following. 
We consider (4.35) and we first perform the summation over /. The quantization condition is 

Qi (9) = mL sinh 8 = 2nl Pl = Q[ (4.36) 

with leN. The converting the sum over I to contour integrals we get 

v v J" d9 F%{8 + in, 9[, flQiffi^fl + in, 8> n , ...,8[) 

jZXlb 2 * Pn(0[,...,9' n ) 

X Q^ mx ( s ' nn 9—mx y\ sinh 9j)—m(R—t) cosh 9— vat yV cosh Oj 1 r^y\ 
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where the contour Cj surrounds the solution of (4.36). When opening the contour to surround 
the real axis in 9, we get the following contribution in the L — > oo limit: 



i r de rde[ de' nj?0l 



F^(9 + m,9[,..., 9 n )F% 1 (9 + i7r,9' n ,...,9' 1 



n! Jc 2vr J 2vr 2vr 



x e 



imx(sinh 9— ■ sinh 9\ Ax— m(R— t) cosh 9—mty^- cosh 9' 



(4.38) 



where the contour C+ is defined as 



d9 
2^ 



(4.39) 



However, there are additional poles of the integrand for 9 = 9'a for j = 1 . . . n, whose contribution 
must be subtracted. 

First we calculate the residue at 9 — > 9[. The behaviour of the form factors is given by the 
kinematical residue equation (2.10). Let us introduce the (partially) connected part of the form 
factor as 



f 



+ m 



i-n^-^) 

3=2 



9-9[ 



(4.40) 



The form factor F sc defined above is only "partially" connected since only one of the singularities 
is subtracted and so it still has poles at 9[ = 9j for j = 2 . . . n. In fact, the singularity structure 
of the original form factor near 9 = 9[ = 9[ is given by 

" / 1 1 \ 

F%(9 + i- ! r,9' 1 ,9[,...)=i(l + J[S(9[ - 9 k ))F°\(02, 9' n ) {^-^ - j-^j (4.41) 

and after subtracting the first pole there remains the second one leading to 

n 

9 2 , ... , 9' n ) = + J] S(9[ - fc ))*?ii(*2, ...,9> n )-— 



+ ... 



k=3 



The connected part satisfies the exchange equation 



F, 



Oi 



?k ■ ■ • > O - ~ 9'k) F n+i,sc(0i\6i' ■■■i^'k^j ■■■■> 9'n) 



n+l,scV"lrl' ■ ■ ■ J "ji 



(4.42) 



(4.43) 



In the case of the crossed form factor one has 



F^ 1 {9 + m 1 9' n ,---A) 



3=2 1 



F 



2 rat 



l[S(9 j -9' 1 ))F^ 1>sc (9[\9 / n ,...,9 2 ) + ... 

3=2 



The residue of the integrand at 9 = 9[ is then expressed as 

g-imi E"=2 sinh 6 j—mR cosh 9^— mt J2"=2 cosh6^. .f^Oi fQ' Q< \^pOi fgt qi •j 

n 

x [( n s ( e i ~ °^ ~ x ) ( imx c ° sh ^ ~ m ( R ~ sinh ^i) + imL c ° sh ^ 



(4.44) 



3=2 



+ iF n+l,c(Q'l\0li ■ ■ ■ >Qn) F n-l(0n> ■ ■ ■ + iF n+l,c(^lWw ■ ■ ■ ^l) F n-li^ ■ ■ ■ >0} 
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There are similar contributions at 9 = 9j for some j > 2, with the role of 9j and 9[ exchanged. 
After integrating over all the 9j one could make a change of variables to obtain the same contri- 
bution n times. However, one has to keep these residues separately because of the poles of the 
connected form factors. Making the change of variables only in the regular terms one obtains 



-imi J]?-2 sinh fl' -mfi cosh 6! -mi cosh 9' p Oi /a' a' \ /a' a' \ 

n 

x [( II S ( e 'i ~ Q 'r>~ X ) ( imx cosh 9 'i - m ( R - l ) sinh °'i> + imL cosh 9 i 
i=2 

The 0{L) term of (4.45) term can be transformed in the L — > oo limit into 

— L_ mL [ ^1 cosh 9[e- mRcosh9 i 
(n-1)! J 2vr 1 

7 — • • • —Fn-M, ■ ■ • • .,9> 2 )e- imx ^ sinhe 'i- mt ^ cosh( 

J 2ir 2ir 



(4.45) 



fit'- ^ F ^ • • • ' O*£l0C ■ ■ ■ , 4) ( 



(n-1)! 



(4.46) 



mR cosh 6 j 



3=2 

x e — imx sinh ®'i ~ mt ^5=2 cosh o'j 

The subtraction of the last term takes into account the exclusion principle 9[ ^ 9'j for j = 2 . . . n, 
which is already present at the level of quantum numbers. The first term in (4.46) gets exactly 
canceled by Z\CQ <n -\ leaving only the second one which is finite as L — > oo. 

The O(L ) terms of (4.45) are regular, therefore it is allowed to replace the summation over 
the rapidities with the appropriate integral. 

Putting everything together, the net result is 

Din =h.S 1 1 fr • • • it F ^ e + ->^> • • • > + ^ c • • • , oi) 



x e 



c+ 

imx(sinh 8— sinh &'A—m(R—t) cosh 8— rat cosh 8'- 



I 1 f dO'i d9!n -imx J2] =2 sinh e 'j ~ mR cosh e 'i- mt Yl]=2 cosh e 'j ( TT qiq' —Q'\ — \\ 
(n- 1)1 J 2vr "' 2vr ^11 1 1 i> > 

x {mx cosh 9[ + im(R - t) sinh 9'JF^ (9' 2 , . . . , O^f-i (C • • • , «£) 

1 /" I r -irnxJ2" =2 sinh ^ - mi? cosh e l- mt E°=2 cosh e 'j 

n\J 2tt " ' 2tt H 

X { F n+l,sc(0'l\0'li ■ ■ ■ i9'n) F n-l(9'ni • • • > #2) + F n+l,sc^'l\^ ■ ■ ■ >^l) F n-l(^ ■ ■ ■ > O) 

9[ O e'j for j = 2..n 



+ 



/§...^i£i(^...,oiei(^...^ 2 )(E e 



(n-1)! 



J'=2 



(4.47) 



4.4 Life is not that simple: L> 2 2 

Now we turn to the evaluation of the 2-particle - 2-particle contribution to the thermal cor- 
relator. The novel feature of this contribution is that the diagonal terms must be separated from 
the non-diagonal ones, since the four-particle form factors (in contrast to the two-particle one 
that appears in Dn) have nonzero residues for the kinematical poles according to (2.10). When 
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evaluated at the diagonal, these singularities are eliminated but result in an ambiguity of the 
diagonal matrix element, which was discussed in much detail in [45]. Once this complication is 
attended to, the evaluation proceeds similarly to that of D12. 

According to the general formalism outlined in section 2.4 we can write 

D 22 = hm D 22 (4.48) 

where 

D22 = C22 - ZxC n + (Z 2 - Z 2 )C 00 (4.49) 

Using 

D n = C n - Z x C m (4.50) 

gives 

D22 = C22 - ZxD n - Z 2 C 00 (4.51) 

The new contribution is 

C 22 = EE ({h,h}\0 1 \{J 1 ,J2}) L ({Jl,J2}\02\{h,l2}) L 

Il>l2 Jl>J2 

x K^\e 1 ,e 2 ;e[,e' 2 ) (4.52) 

where 

imx(sinh 0i+sinh 82 — sinh 8[ — sinh 2 )p — m(R— t)(cosh #i+cosh 82) 
t,x v 1) 2i V\i "2) e " 

X g— mt(cosh0j+coah02) (4 53) 

We can separate the sum into diagonal and non-diagonal part: 

E E = E( termswith { J i' J 2} = i / i' / 2})+ E E ' ( 4 - 54 ) 

h>hJi>J2 h>h h>l2.Ji>J2 
where the prime means that {Ji, J2} 7^ ^2}- 

4.4.1 Evaluating Z 2 

First of all, we need the two-point contribution to the partition function. This is easy to 
obtain: 



^ 2 _ e -m-R(cosh6»i+cosh6» 2 ) J ^ J e -mi?(coshe 1 +cosh0 2 ) (4 55) 

(4.56) 



2 ^— ' 2 



To convert the sums to integrals, we need the two-particle density of states 

P2(0i,O 2 ) = m 2 L 2 cosh 0i cosh Q 2 + mL(cosh0i + cosh 2 M#i - 2 ) (4.57) 

and also the density of states on the diagonal I\ = I 2 which can be obtained as the derivative of 
the degenerate (0i = 02) Bethe-Yang quantization condition as follows 2 

Qd(8i) = mL sinh 0i = 2-7r/i 

Q' d {6i) = mLcosh0i = pi(0i) (4.58) 
i.e. it coincides with the one-particle density p\. The result is 

Z 2 = \ [^ p2{ 8 l ,62)e- mR ( cM+cosh ^-l ff- Pl (9)e- 2mRcoshe (4.59) 

Z J Z7T J Z7T Z J Z7T 

2 Note that due to (2.19) 5(0) = 0. 
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i [ dOi f de 2 e - mR ( cosh ^+ coshe ^ 



4.4.2 The diagonal sum 

Using (2.26), the diagonal matrix element is 

{{h,h\\0\{h,h}) L = (a a \ ( 4 - 60 ) 

Substituting these into the diagonal sum 

({h,h} {h,h}) L ({Iuh} \0 2 \ { Iul2 }) Le -mR(coshe 1+ co Sh e 2 ) (4 M) 

h>h 

we have terms that can be ordered by the number of F± s factors they contain. The term which 
contains two F^ s factors can be written as 

^ g— mii(cosh 0i+cosh 82) 

E n (ft \2 F 4s 1 (9x,9 2 )F 4s 2 (9 1 ,9 2 ) = 

— mi?(cosh 61 +cosh 82 ) 

— —jh~a~v 2 — F^ s i (9 1 ,9 2 )F^ s 2 (e 1 ,e 2 ) = 

82) 

2nJ 2tt p 2 {9 u 9 2 ) — F ^^ 2 )F^{6 u e 2 ) 
, d g i e - 2mfl (cos hei ) mLcosh g 
~J 27 p~2~(hM 2 43 {ei >° l)F ^ { ° M (4 ' 62) 

where we added and subtracted the diagonal 9\ = 6 2 . We used that the density of diagonal 
(#1 — #2) two-particle states is given by mL cosh 9\. The first term is 0(L~ 2 ) while the second is 
0(L~ 3 ) and so they vanish as L oo. 

There are two terms containing a single F± s . One of them is 

1 — mi?(cosh 0i+cosh 82) , > 

9 E uTTxi F^(6 1 ,9 2 ) (p 1 (e 1 )F^+ Pl (e 2 )F^+ P2 (e 1 ,e 2 ) (0 2 )\ (4.63) 

and the other can be obtained by interchanging 0\ and 2 . In writing the above formula we 
already included the I\ = I 2 diagonal, using again that it is suppressed by an L" 1 factor. For 
L — > oo we get 

1 r d^dfo e _ mR{M + msiie2)F o m (4 g4) 

2 J Z7T Z7T 

Similarly, the other term yields 

I J ^^ e ~ m R(cos h 8 1+ cos^82) F O2 {0ue2){Oi) (4 g5) 



The terms without F^ s give 

J?,(cosh 

pWM 



— mi?(cosh 0i+cosh 82) , 

E " — n( o w\2 — pi(^i)^ 1 +pi{02)fZ 1 +P2(e 1 ,e 2 ){o 1 ) 



h>h 



x (pi(#i)F£ 2 +pi(9 2 )F£ + p 2 (9 1 ,9 2 ) (0 2 )) (4.66) 

We can replace 

E - \ E -| E 

Il>il2 /l,/2 Jl=-f2 
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and after converting the sums to integrals we obtain 

1 ^e--.-™, tpm ^ +pm ^ +p2{eue2){0l 



2 J 2ir J 2tt p 2 (9i,e 2 ) 
Pi(0i)f£ 2 + pi(9 2 )F^ + p 2 (e 1 ,e 2 ) (C 2 ) ) - diagonal term 



fde 1 


fd9 2 


1 2vr J 


2n 



1 -mR(cosh 0i +cosh 2 ) pOi 77O2 (cosll 01 + COsll 2 ) 2 



-e 



2" -2c 2c cos h 6»i cosh l9 2 

.01 



+ mL cosh eie- mR(coshei+cosh02) ((d) F° 2 + (O a ) F 2c 

+ I e -mfl(cosh0 1+ cosh0 2 ) ( m 2 L 2 cosh0iCosh ^ + mL(cQsh ^ +cosh ^ 2 ^^ i _ ^ ^ ^ 

\ J ^~ 2mRcosh9 [mLcoshe (d) (0 2 ) + 2 ((d) + (0 2 ) (4.68) 

where we dropped terms that vanish as L — > oo. This has terms which diverge in the limit; 
however, we must now add the "counter terms" 

-Z 2 C 00 = - f— mLcoshee~ 2mRcoshe {d) (0 2 ) (4.69) 
2 J 2ir 

2 \vr?I? cosh Q\ cosh 9 2 



and 



2 J 2vr J 2ir 

+mL(cosh0i + cosh 2 M0i - e 2 )]e- mR{coshei+coshd2) (C»i) (0 2 ) 



-Ztf^f ^i F ^(e 1 + i 7r,e 2 )F^(9 1 ,9 2 ^ 

^ gimx(sinh 0i— sinh 02 )g— m(R— t) cosh 0i ^— mi cosh 02 

-Z a y ^ (i^ 1 (0 2 ) + F^ 2 (Oi)) e" mRcoshe (4.70) 
These cancel all the divergences leaving us with the final result for the diagonal contribution: 
D (diag) = 1 J ^^ e --^^+ cosh ^)(FfH^i^2)(d) + F^(0 1 ,0 2 )(d)) 

1 [ dfi r dd 2 mfl ( cosh 1+ cosh02) F 0i pQ 2 (cosh 0i + cosh0 2 ) 2 

2 ] 2W 2i 2c 2c cosh 0! cosh 2 

- J ^- e ~ 2mRcoshe ((d) F^ + (0 2 ) F^ 1 ) (4.71) 
However, the first piece of —Z\D\\: 

/ x,V ft„\-p0 2 (a a , • \ Jmi (sinh 0i -sinh 2 ) -m(R-t) cosh 0i -mt cosh 2 



-Zx — — F 2 ^(0i + i7r,9 2 )F; Y 2 (0i,02 + i 7r ) e lma: ^ innWl - slimt '^e- m ^- t J cosnei e 

7 27T 7 27T 

(4.72) 

is not canceled by the diagonal part. We now turn to the evaluation of the non-diagonal contri- 
bution, which does eliminate this last divergence, as explicitly demonstrated in appendix C.2. 

4.4.3 Evaluating the non-diagonal part 

Now we must evaluate 

^(nondiag) = ^ ^ ' , / 2 } | ^ | { Jl , J 2 }) L ({ Jl , J 2 } | d | Ul , / 2 }) L 

h>h Ji>J2 

xKQieufaOlX) (4-73) 
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where 

O ■ O' O' \ ima;(sinh 6*i+sinh 82 — sinh 9^ — sinh 8' 2 ) — m(R— t)(cosh 6*i+cosh 82) 

n t,i lyl) "2t u \i u 2) ~ e e 

^g-mijcoshej+cosh^) (4 74) 

Using (2.24), the matrix elements are of the form 



({I 1 ,L 2 }\0 1 \{J 1 ,J 2 }) L 



F^ 1 (6 2 + iTr,e 1 +iTr,e[,9' 2 ) 



2 



{{JuJ 2 }\0 2 \{h,I 2 }) L - p2{e ^ 2Y/2p2{e ,^ 2Y/2 (4-75) 
The quantization conditions read 

Qi(6>i,0 2 ) = mL sinh 9 1 + 5{9 1 - 9 2 ) = 2nh 

Q2(0i,e 2 ) = mLsmh9 2 + 5(9 2 - 9 1 ) = 2ttI 2 (4.76) 

and 

Q v (9' 1 ,9' 2 ) = mL sinh 9[ + 5(9[ - 9' 2 ) = 2vrJi 

Q 2/ (9[,9' 2 ) = mL sinh 0£ + 8{9' 2 - 9[) = 2ttJ 2 (4.77) 

Now we can write 

({I 1 J 2 }\0 1 \{J 1 ,J 2 }) L {{J 1 ,J 2 }\0 2 \{L 1 J 2 }) L K^\e 1 ,9 2 -,9' 1 ,e' 2 ) (4.78) 

/[' d9' d9' 
f ~2~t 2**** { ° 2 + i7F ' ° l + ™> 9 '^ F ° 2 ^ + i7r ' ^ + i7r ' 

( e iQi'(fi>^) + 1) (e^aWa) + l) k ^ 

where Cj x j 2 is a multi-contour surrounding the solution of (4.77). When we open the multi- 
contours to surround the real axes we encounter new singularities. These can be classified as 
follows: 

1. QF-singularities: partly from the Qs, partly from the Fs: 

9[ = 9 1 and Q 2 ,(9' 1 ,9' 2 ) = 2ttJ 2 (4.80) 

9[ = 9 2 and Q 2 ,(9[, 9' 2 ) = 2tt J 2 (4.81) 

9' 2 = 9 1 and Qi>(9[, 9 2 ) = 2tt Ji (4.82) 

9' 2 = 9 2 and Qv{#y, 6 2 ) = 2ir Ji (4.83) 



2. i^F-singularities: come from the Fs 



A — 9' 2 — 9\ 
9'i = 9' 2 = 9 2 



(4.84) 
(4.85) 



3. (Spurious) QQ-singularities: these result from 



and 9' 2 = 6 2 (4.86) 
' - ° and 9' 2 = 9 1 (4.87) 



V X = (71 

9'i = 9 2 

(it turns out that eventually these do not give any contributions in the L — > oo limit). 

Since the evaluation of these contributions is the same as for D\ 2 , the details are relegated to 
appendix C. The upshot is that 

^,11 d9\d9' Iff M d9' 2 .„„ , , ,„ . 

where C is the open multi-contour. Here we also used the fact that there are no singularities at 
J\ = J 2 because the form factors vanish. 
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4.4.4 End result for D 22 

Putting together the results (C.15), (C.22) and (C.29) of appendix C with (4.71) one obtains 

(4.89) 



D22 
1 
4 



d9id9 2 ff d9' l d9' 2 . , a ,-, I? o 2l 



+ 



+ 



2vr 2vr J J c++ 2vr 2vr 

\ J [ ^f*e^^*+«**> (f£(9 u 9 2 ) (0 2 ) + ffi-fa, (9 2 ) (Ox> 

1 /"/■ dgi dgg m i?(cosh ex +cosh0 2 ) F Oi F 2 (cosh^i + cosh0 2 ) 2 

2 77 2vr 2vr 2c 2c cosh X cosh 2 



2 + iir,9' l + m,9i,9 2 ) 



xc 



imx(sinh 02- sinh Sj)^,— mi? cosh $i —m(R—t) cosh #2 „ — TOt cosh 2 



((mxcosh^i - imi sinh 0x)(l - S(0 2 - 6 1 )S{9 1 - 9 2 )) + y>(0i - 9' 2 )S{9'2 - 9 1 )S{9 1 - 9 2 )) 



+ 



d9 x d9 2 f d9' 2 



2ir 2ir 



2tt 



3 imx(sinh 82 — sinh 8' 2 ) p —mR. cosh 8\ —m(R—t) cosh 82 — mt cosh 8' 2 



e v 2 'e x e v ' ^e 

x (F 2 Ol (0 2 + wr, 9' 2 )Fgl{9' 2 , 9 1 \9 1 , 9 2 ) + F°*{9' 2 + in, 9 2 )Fg\{9' 2 , 9 X \9 X , 9 2 ) 
+ (0i O 2 ) 



gimo;(sinh 6*2 — sinh 61 )^— m(i?— t) (cosh 6*i+cosh 82) ^—2mt cosh 0i 



02)F 2 Ol (9 2 + i7r,9 1 )F 2 ° 2 {9 1 + iir,9 2 ) 





d9 2 - 


HZ 


2vr . 



x5(0i 

+ (01 O 2 



d0id92 FOl 
2ir 2vr 2 



F 2 ° 1 (9 1 + m,9 2 )F 2 ° 2 {92 + i7r, ^ x ) e ima; ( sinh i ~ sinh <? 2 )-^(2i?-t) cosh 6 1 -mt cosh 



02 



where the function is defined in (C.16), and C++ denotes the integration contour specified 
in (3.31). All the other integrals are taken over real values of their variables. 

4.5 Discussion of the proposal of LeClair and Mussardo 

In [34] LeClair and Mussardo introduced a regularization scheme for finite temperature corre- 
lation functions. The two main assumptions of the proposal are that the spectral expansion should 
be built using the zero-temperature form factors, and that the only effect of finite temperature 
is an appropriate modification (dressing) of the statistical weight functions and the one-particle 
energies and momenta. In the case of one-point function the proposed formula was proven to 
be correct up to the third order in [45]; an all-orders proof is also possible [46]. However, the 
two-point function seems to be more problematic [35, 43]. 

In the following we compare our results to the proposal of [34]. For the two-point functions 
their formula reads 



<0(x,i)0(O,O)) 



R 



«o> R ) 2 + E^ £ / 

N=l ' Oi=±l J 

x|(O|0|0 1 ...0 iV > CTi ... CTjv | 2 



d9i d9 N 
~2~k " ' ~2tT 



N 



(4.90) 



where f a 



1/(1 + e~ a ^ e ^), £j = e(9j)/R and kj = k(9j) with e{9) being the solution of the 
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TBA equation 

dff 



r do 

e{6) = mR cosh 9- / — tp(9 - 9') log(l + e~ 6 ^ ) (4.91) 
and /c(#) is given by 

k{9) = msmh9 + J d9'5(9 -9') Pl (9') 
2^i(0)(l + e £(e) ) = mcosh0 + I d9'tp{9 - 9') Pl (9') (4.92) 



The form factors appearing in (4.90) are defined by 

(O|0|0i . . . 9 N ) ai ^ N = F^(9 1 - iTrorj, . . . , 9 N - ma N ) Gj = (1 - ^)/2 e {0, 1} 

The interpretation of the series (4.90) is as follows: the excitations with Uj = +1 or <jj = — 1 
represent particles or holes over an infinite volume thermal state. Therefore the statistical weight 
functions and the one-particle energies and momenta are given by the dressed values as calculated 
in TBA. 

If the (4.90) series expression were correct, then a systematic double expansion in terms of 
e —mt anc j e -m(-R-t) should reproduce our results. Indeed, the first few terms indicate that this 
might be true. It was already pointed out in [57] that to the lowest order in e~ mR the term N = 2 
reproduces our Z?o2 + -Cn + -^20- Moreover, in the case of our D\2, the last line of formula (4.22) 
suggests the dressing 

p —mt cosh 9+imx sinh 6 

g— rat cosh 9+imx sinh 9 y _ , Qi^—2mR\ 

_|_ g— mR cosh 9 ' 

for the exponential factors in Z?oi- This pattern repeats itself and similar contributions can be 
found in D\ n (4.47), which suggest a dressing of the factors in Do,n-i- However, the situation is 
more complicated as we consider higher order terms. 

First of all observe, that the formula (4.90) is not well-defined for > 3. There appear 
second order poles whenever the rapidity of a particle approaches the rapidity of a hole, and in 
the original work [34] it is not explained how to integrate over these singularities. Note that it 
was the evaluation of these ill-defined terms which required a lot of effort in our evaluation of the 
two-point function. 

Based on the form of our results (4.22), (4.47) and (4.89) it seems unlikely, that any regularized 
form of (4.90) would be correct. However, at present we cannot make any definitive statement 
about this issue. The inspection of higher order terms (D nm with n, m > 2) might decide whether 
there exists a neat formula for the two-point function, possibly with a structure similar to (4.90) 
but with different dressing prescriptions. This problem is left for future work. 



5 Second order form factor perturbation theory 

As a further application of the framework presented here, we show how to derive the main 
results of the paper [50] on second order form factor perturbation theory using the present for- 
malism. We simplify the presentation by considering a theory with a single massive particle in 
its spectrum instead of the double sine-Gordon theory treated in [50]. Consider modifying the 
Hamiltonian of an integrable model as follows: 

^nonintcgrablc -^integrable ^ J dx^(t, x) (^'1) 

where ^ denotes a local (Lorentz scalar) field which breaks integrability. Corrections that are 
first order in A were derived in [51], but when evaluating the second order one encounters the 
same difficulties with disconnected terms as in the case of the thermal two-point function. The 
principle of the solution to this problem is the same as for the thermal correlator: we perform 
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perturbation theory in finite volume, express the quantities we are interested in and then take the 
limit L —7- oo. In the approach of [50] it was necessary to compute some part of the discrete sum 
over the finite volume quantum numbers explicitly; we show that this can be greatly simplified by 
applying the residue methods of the present work. 

The general perturbation theory formula for second order corrections to energy levels is 

'{i\Hi\k)\ 2 



5E, = £ 



E 



(o) 



E 



(o) 



, H 1 = A J dxV(t,x) 



(5.2) 



kj^i i ^k 

therefore the correction to the vacuum level can be written as 



SE Q 



k^O 



: expif (p(0, 0) : 



k 



E (0) _ E (o) 



(5.3) 



J k ^0 

The summation goes over all excited states in the spectrum (with zero total momentum selected 
for by translational invariance), which can be described using the Bethe-Yang picture of section 
2.3. The leading contribution is given by the state containing a single stationary particle, and can 
be written as 



X 2 L 2 



SE (L) -- 

Using the relation (2.24) we obtain 



: exp if 99(0,0) : {0} 
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+ O (e"^ L ) 



5E (L) 



-X 2 L 2 



*l 2 
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#|2 



O (e-» L ) = -X 2 L ] —^- + O (e-^ L ) 



p%(0)m 

which results in the following shift of the bulk energy density 
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-A 



9\2 



(5.4) 



(5.5) 



(5.6) 



Next we are interested in the correction to the particle mass. This can be obtained by evaluating 
the correction to the first zero-momentum excited level in the finite volume system and then 
taking the limit 

(5.7) 



6m = lim 5E 1 (L) - 6E (L) 

L— >oo 



The correction to the first level can be written as 
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F 3 (iir,9,-9)F 3 (9, 



, in 



pi(0)p 2 (0, -9) (2m cosh 9 - m) 
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(where we omitted the states with three or more particles) where 9 is the solution of (cf. subsection 
2.3). 

Q{9) = mLsinh 9 + 5(29) = 2irl , IgN+J (5.9) 



and 

p 2 (0, -9) = mL cosh 9(mL cosh 9 + 2</j(20)) 
Extending the sum over 9 to negative values and performing the residue trick we get 



(5.10) 
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Cl 2tt e i( 2W + 1 



m 3 (2 cosh 9 — 1) cosh 9 



(5-11) 
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where 

p 2 {9) = mL cosh 9 + 2^(20) (5.12) 

is nothing else than the density of two-particle states with zero total momentum. Using (2.10), 
the form factor has the following singularity at 6 = 



|F 3 (i7T,0,-0)| ! 



16[Fi| 



+ O(0°) 



(5.13) 



where we also used 5(0) = —1. Subtracting and adding the singular term at the origin results in 
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16|Fi| 
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27r e^C) + 1 m 3 sinh 2 9 cosh 



(5.14) 



We can then open the contours to surround the real axis. In the first term, only the upper contour 
contributes in the infinite volume limit, and one can also take e — > 0. In the second term, we 
can open the contour, but we must also subtract the contribution of the double pole at the origin 
since that is not included in the original sum: 
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(5.15) 
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The first integral is 



de 



2tt e*W) + 1 m 3 sinh 2 9 cosh 9 



16|Fi| 



A 2 



27T m 3 sinh (9 + ie) cosh(0 + ie) 



de 



1 



16|Fi| 



Co 



2vr e*0( e ) + 1 m 3 sinh 2 (9 cosh (9 



A^ 
2 



d0 



16 |Fi 
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while the second integral is given by the residue theorem as 
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Finally, using (5.7) and (5.5) we obtain the mass correction 
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AM f 



iF 3 (m,0,-0)r 



16IF1 



2n I m 3 (2 cosh — 1) cosh m 3 sinh 2 cosh 9 



+ A' 



16|Fi| 

4m 3 



[l-<p(0)) (5.20) 



which agrees with the result in [50] . Note that the leading bulk term drops out from the difference 
of the energy levels as it indeed should. 

The contributions of higher-particle states to the spectral sum can be computed analogously. 
It is straightforward to verify that the n particle state term always contains a double pole part 
analogous to the one treated above, which exactly cancels the n — 1 particle contribution to the 
vacuum level in (5.3). 
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6 Conclusions and outlook 



First let us sum up what has been achieved in this paper. Using the idea of finite volume 
regularization and multi-dimensional residue techniques we have developed a systematic technique 
to evaluate the form factor expansion for the finite-temperature two-point function in integrable 
field theories. In fact, as the examples of the zero-temperature three-point function and of form 
factor perturbation theory show, the approach can be applied to many problems involving spectral 
sums with singularities coming from the presence of disconnected terms. Albeit it was expected 
on general physical grounds, it is an important fact that our calculation demonstrated that the 
resulting expressions for the correlators (and for the mass gap in the case of FFPT) are well-defined 
when removing the regulator by taking the infinite volume limit. 

For the three-point function it is apparent that the resulting formula is just the proper way of 
separating the disconnected pieces and can in fact be written down directly by inspection of the 
infinite volume expression. This is due to the disconnected pieces appearing linearly. However, 
the other two cases involve the disconnected terms squared, similarly to the case of one-point 
functions of bulk operators with boundaries [49]. While it is well-known that the resulting terms, 
naively containing squares of Dirac 5 functions can be regularized in a finite box, the correct result 
can only be obtained by carefully taking into account that the finite volume spectrum is different 
from that of a non-interacting system. This was pointed out for the one-point functions in our 
previous paper [45], and is manifested by the explicit dependence of the results (4.22), (4.89) and 
(5.20) on the S matrix (directly or via the derivative (p of the phase-shift). The LeClair-Mussardo 
proposal for the one and two-point functions tries to capture this feature by a TBA dressing of the 
energy and momentum of the finite-temperature quasi-particles. Contrary to the one-point case 
[45], our result for the two-point function does not confirm their conjectured expression (which 
is, in any case, eventually ill-defined). Despite some partial indications of resummation, it is not 
obvious whether the interaction dependence can be summed up to yield some simple dressing 
prescription. 

This leads us to one of the main open questions, namely, to investigate the possibility of such 
resummation and find out whether there is a way to introduce some sort of dressing prescription 
to simplify the expansion by systematically combining contributions. At this point this seems to 
require the evaluation of higher orders, which is in principle straightforward, but an extremely 
tedious task. Therefore another important (albeit technical) issue is to simplify the method of 
evaluating the contributions to the expansion. 

It is also very important, especially in view of potential applications, to extend the method to 
non-diagonal theories. While this is in principle straightforward (for the basic ideas cf. [53] in the 
framework of 0(3) model), it would be desirable to have some efficient approach to characterizing 
the finite-volume form factors of non-diagonal models for general number of particles. Work in 
this direction is in progress. This is even more important, since at present the most we can show 
for testing the expansion for the thermal correlator is its internal consistency. Consistency is 
shown by two facts: (1) that terms divergent for large volumes cancel in the final result order- 
by-order and (2) that different orders of performing the summation lead to identical results, as 
demonstrated in for Dyi- These are indeed very nontrivial tests of the calculation, but a physical 
application of the method would be much better. 

Another important issue, especially in view of potential applications to non-relativistic sys- 
tems along the lines of [58], is the extension to include a nonzero chemical potential. At present 
it is not entirely clear how to do that, but rewriting the expansion through some partial resum- 
mation/dressing procedure could be helpful (in analogy to the way the dependence is introduced 
into the LeClair-Mussardo formula for the one-point function [34]). 

Finally we comment on the relation of our results to the recent work by Essler and Konik 
[53]. Their finite- volume calculation is essentially the evaluation of Dyi using the one-particle 
summation, with an explicit summation of the discrete part instead of a residue trick. However, 
this approach is very limited: it can only be applied to D\ n since all other contributions require 
summation over states with two (or more) particles and there is no obvious way to perform the 
discrete sums directly. This is where the multi-dimensional residue method presented here is 
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so powerful since it makes the evaluation of such sums a mechanical (albeit somewhat tedious) 
exercise. 

The infinite- volume regularization method of [53], on the other hand, is plagued by (at least 
potential, but most likely actual) ambiguities: for states containing more than one particle, regu- 
larization by point-splitting in rapidity space is ambiguous (direction-dependent) at the locations 
in rapidity space where the form factors have either of the two types of disconnected contributions 
described in subsection 2.3. These ambiguities were analyzed in much detail in our previous paper 
[45]. Therefore (at least at the present state of art) the only safe method to evaluate multi-particle 
contributions is by finite volume regularization, and the only systematic way to perform the sum- 
mations is by using the multi-dimensional residue theorem, i.e. in the framework presented here. 
As mentioned above, however, it is an important goal to simplify the method of calculation, which 
could potentially lead to dispensing with these technical requirements in the end. 
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A Multi-dimensional residue formula 

Suppose that we have functions g(z), fi(z), . . . , f n (z) of n complex variables z = (z±, . . . , z n ). 
Let us take a multi-contour C in C n (i.e. a direct sum of elementary multi-contours that are 
defined as direct products of n one-dimensional contours; without loss of generality - due to the 
linearity of integration - we may suppose that it is a single (i.e. monomial) product contour 
C = C\ X ■ • • X C n ). Let us suppose that the equations 

fk(z)=0 k = l,...,n (A.l) 

have a single solution z* — (z*i, . . . , z vn ) such that for each k, Zj~* is inside C\~. Then we have the 
formula 

dzi dz n g(z) g(z*) ^ 



c 2m"' 2m h{z) . . . f n {z) d t ( dp, 

provided that the determinant does not vanish (this will always be the case in our calculations). 
Note that if any of the fj. is nonzero everywhere inside its contour Cfe then the integral vanishes 
since Ck can be shrunk to a point. Therefore, all the usual contour deformation arguments work 
as long as the contour deformations take place away from the analytic variety defined by 

de«(f)=0 (A.3) 
B Finite volume FF and phase conventions 

One has to evaluate the products of finite volume form factors. It follows from the crossing 
formula (2.5) that 

{{h . . . I N }\Oi(0)\{Ji . . . Jm}) l ({Ji ■ ■ ■ J M }\O 2 (0)\{h ■ ■ ■ In}) l = 

ggWgi + i*,...,8 N + m,e' M ,..., 0'i)Fn+m( S i + ; ; ; ; 9 'm + ^ d N, ■ ■ ■ , gi) (B.l) 

PN (0 1 ,...,e N )pM(0' 1 ,..-,o' M ) 

In the following we show that in unitary models the crossing procedure described by (B.l) 
reproduces the usual complex conjugation of the matrix element. In unitary models the phase of 
the form factor is given by (allow for an extra sign ambiguity) 



F N (e 1 ,...,e N )= F N (e 1 ,...,e N ) x JlsiOi-ej) e«eM (b.2) 



i<j 
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In particular F\ is always real. The extension to include also "bra" vectors reads 
F N+M (0[ + iir, . . . , 9' M + i-K, 9i, . . . , 9 N ) = 



F N+M (e' 1 + iTr,...,e' M + m,e 1 ,...,e N ) x /n^-^) * JH s ( 0, k - d i) (b.3) 

V i<j V k<l 



6i 6 K 

The complex conjugation property is then given by 



F N+M (9[ +iT,...,6 , M + iir,0i,..., 9 N )) = F N+M (9' M + ix, . . . , 0[ + itt, n , . . . , 9{) (B.4) 



This way one can avoid the operation of complex conjugation and one can work with analytic 
functions in the complex plane. 

C Evaluating the subtractions in the non-diagonal part of D22 

We first substitute 



E-lE-lE (0.1: 



2 ^ 2 

C.l Spurious QQ singularities 

Let us consider the first family of such singularities, which is when 

6[ = 9 X and 9' 2 = 9 2 (C.2) 

The evaluation of such terms is complicated by the fact that the diagonal limit of F4 is undefined 
(i.e. direction dependent). Fortunately, after converting the sums to integrals a factor of I//J2 
remains therefore all such terms are of order 0(L~ 2 ). Similar considerations apply to the other 

case 

9[ = 9 2 and 9' 2 = 9 X (C.3) 

C.2 QF singularities 

We consider the case 

0i = 0i and Q 2 '(9' 1 ,9' 2 ) = 2vrJ 2 (C.4) 
Just as in the case of D\ 2 , the contribution can be split into a double and a single pole part. 

C.2.1 The double pole part 

The double pole part is given by 

n v 1 f d9^f d9[ Kff(9 1 ,9 2 -9' 1 ,e' 2 ) 

QF P2{0i,e 2 ) J Cj2ll 2vr f 01 2vr (jQAe'M + 1) (jQAW + l) 



1 (1 _ S {9 2 - 9 1 )S(9' 1 - 9' 2 ))F°i(9 2 + m, 9' 2 ) 



n - &i 

i 



(1 - S(9' 2 - 9[)S(0 1 - 2 ))F 2 ° 2 (0 2 + in, 9 2 ) (C.5) 



We have 

On the other hand 

e imLsinh6l S{9 1 -9 2 ) = 1 (C.7) 
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l.C 

c 



iQAh,*,) = s(e l - e' 2 )s(e 2 - Ox) (c.j 

Using S(9)S(—9) = 1 we can rearrange the contribution as 



/l>/2 



x " 



i ( ( i _ S (d 2 - e 1 )s(e[ - e' 2 )) + (i - s(9' 2 - e' l )s{e 1 - e 2 ))) (c.9) 



We now use 

and after manipulations similar to those in subsection 4.2.2 we arrive at 

»qf = -E-^rxyi ^to + ivA^M + i*,**) (en) 

x / n ffl ^ n ( ( mx cosh 01 " sinh ^X 1 - 5 ( 2 - - 6 2 )) 



( e iQ 2 ,(0i,0 2 ) + ij 
+ip(6 1 - 9' 2 )S{9' 2 - 9 1 )S(9 1 - 9 2 ) + mLcosh^i) 

The full contribution is then obtained by adding the three other double pole terms pertaining to 
the other QF-singularities in (4.83), which can also be obtained by suitably permuting the rapidity 
variables. Converting the sum over 6\, 9 2 to integrals, taking care to subtract the diagonal 9\ = 9 2 
we obtain for the full double pole contribution the expression: 

Dqf = fJ^f^F 2 °^9 2 + l ^)F°i(9> 2 + l n,9 2 ) 

imx(s'mh92— sinhSj)^,— mRcoshBi — m(R— t) cosh 82 — mt cosh 8' 2 

x ((mx cosh 9 X - wnisinh0i)(l - S{9 2 - 9 l )S(9 1 - 9 2 )) 



+y{9i - 9' 2 )S(9' 2 - 9 1 )S(9 1 - 6 2 ) + mL cosh 9i 
Jwl ^ F 2° 1 ^+i^2)F 2 ° 2 (e' 2 + in,9 1 ) 



^ gimx'(sinh 9\ —sinh 6' 2 ) ^—m(2R—t) cosh B\ g~ m * cosh 9' 2 -^2) 

Now we recall the leftover counter term from (4.72) 
-Z x j ^ J ^F 2 Ol (9 1 + i7r,9 2 )F 2 ° 2 (9 u 9 2 + in)e imx{sinh01 ~ sinhd2) ^ 

(C.13) 

with 



mL J ^cosh9e- mRcoshe (C.14) 



35 



and see that it exactly cancels the 0{L) part, leaving us with the finite expression 
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2 + Z7T,0l, 
mt cosh 6o 



(C.15) 



C.2.2 Single pole contributions 

Once again, we consider the 9[ = 9\ case and introduce the notation: 



F° 1 {e 2 + m,e l + m,9' l ,e' 2 ) 



(1 _ 5 (0 2 - 0i)S(0i - 2 ))if J (0 2 + ivr,0 2 ) 



n -1 



(C.16) 



and similarly 



Ff 2 (0£ + iMi+iMi,0 2 ) 
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■(1 - 5(0 2 - 01)5(0! - 2 ))if 2 (0 2 + tTT, 02 



+ F 4 ^(0 2 ,0' 1 |0 1 ,0 2 ) 



(C.17) 



The contribution has the following form 



ell 
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2f 01 2vr f Cj2h 2vr ( e «Ji'(^) + l) (e 
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(1 _ 5(0 2 - 005(0; - 02))F 2 Ol (0 2 + i7r,0 2 )F 4 ^(0 2 ,0i!0l,0 2 ) 
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Recalling (C.8) and evaluating the residue integrals 
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oil 
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d0 2 
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(C.19) 



Note that this vanishes when 0i = 2 due to the form factors vanishing, so when putting in the 
01 , 02 summation we can include the diagonal 9\ = 9 2 . Converting the summation to integrals we 
obtain 

gll - f f ^1^^. f imx(smh 8 2 -sinh 9' 2 ) -mR cosh 9 1 -m(R-t) cosh 82 -mt cosh 8' 

QF 4 J J 2vr 2vr J 2ir 

x (if 1 (0 2 + iTT, 0' 2 )Fgl(9' 2 , 0i|0i, 2 ) 

+5(0 2 - 0i)5(0i - 2 )if 2 (0 2 + *7r, 92)Fgl(9 2 , 9^, 9' 2 ] 



(C.20) 
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or, using the definition of F^ sc and the form factor equation (2.8) 

gll - [ f f imxjsmh 6 2 -sinh 6i, ) c ~mR cosh 9 1 -m(R-t) cosh 6> 2 c ~mt cosh 6' n 

QF 4 J J 2vr 2vr J 2vr 

x (f 2 Gi (0 2 + iTr, ^)F 4 ^(^,e 1 |(9 1 ,(9 2 ) 

+F 2 ° 2 (e 2 +m,e 2 )Ff&{e , 2 ,e 1 \e 1 ,e 2 )) (c.21) 



5QF 2 / / 2vr 2vr / 2vr 



The full contribution is then obtained by adding the three other single pole terms pertaining to 
the other QF-singularities in (4.83), which can also be obtained by suitably permuting the rapidity 
variables: 

\irnx (sinh 2 -sinh 6' 2 ) e ~mR cosh B x Q -m(R-t) cosh 9 2 ^-mt cosh 6' 2 

x [f°-{9 2 +m,e' 2 )Fgi{9' 2 ,e 1 \GiM + F^(e' 2 + in,e 2 )F^ c (e' 2 ,e 1 \e u e 2 )] 
+ (e 1 ^e 2 )] (c.22) 

C.3 FF singularities 

First let us consider the 8^ = 9' 2 = 9\ case. We need to separate the singular terms from the 
form factors, fro which we introduce a new function F^dc defined by 

F^(8 2 +m,8 l +m,8[,9' 2 ) = -±—(l-S(9 2 -9 1 )S(9' l -9' 2 ))F 2 °^9 2 + iir,9' 2 ) 

til — ti x 

+ r^^^ " ^ ~ S( - 02 ~ ^)) i? 2° 1 (^2 + in,9[) 
tii - ti 2 

+ F°l(9 2 ,9i\9[,9 2 ) (C.23) 

and similarly 

F°*(8 2 +m,8[+m,8 u 9 2 ) = — -i— (1 - S{9' 2 - 9[)S(9i - 9 2 ))F°*{9' 2 + m, 9 2 ) 

f i — #1 

+ W^¥i {s{e ' 2 " e ' l) ~ s{91 " e ^ F ° 2 ^ + m ^ 

+ ^OMl^fc) (C24) 

Only the cross terms can contribute, otherwise at least one of the contour integrals can be shrunk 
to a point. We obtain 

Ql 1 1 / d9' 2 I d0[ K^(9i,9 2 ;9[,9' 2 ) 

FF ~ ~ 2_> 



7 ^ 2 P2(0i,02) 2 J 9l 2vr J 6l 2vr ( e *Ji'(«i,«J) + l) ( e ^ 2 'M) + i) 



1 1 

" '?2 j 



-F 2 Ol ^2+ivr,^)F 2 O2 (0 2 + ivr,0 2 



$i — 8'y 8\ — 9' 2 

x2(l + S(8i - 8 2 ))(l + S{9 2 - 8i)) (C.25) 

We have 

J£^ty (J)^ 9 2 ' 9l 9l) g*™a;'(sinh 02— sinh 0i)g—m(/?—t)(cosh 0i+cosh02)g—2mt cosh 9i (C 26) 

and 

jQAeifa) = gimisinhS! = s(9 2 -9i) (C.27) 
e iQAW = e imLail0hdl = S(9 2 -9i) (C.28) 
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The contribution S FF from 9[ = 6' 2 = #2 can be obtained by interchanging 6\ and 62- Putting in 
the 9\ and 62 integrals we obtain 



Sff 



d(h f d9 2 



gima;(sinh 02 — sinh m(R— t)(cosh (?i+cosh 82) q— 2mt cosh ( 



x S(0i - 2 )F^(9 2 + m^ x )F 2 ° 2 {9 1 + wr, 2 ) + (#1 fc) 



(C.29) 



(in this case the diagonal subtraction is 0(L 1 ), so it does not give a contribution in the infinite 
volume limit). 
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